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I.        INTRODUCTION 

From  July  15th  to  August  17th  of  1997,  trial  runs  using  1:24  scale  models  of  a  T- 
ACS  Auxiliary  Crane  Ship,  a  DDG-963  Class  Ship,  a  Commercial  Container  Ship  and 
Lighter  Barges  were  conducted  in  the  seakeeping  basin  of  the  David  Taylor  Model  Basin 
at  Carderock  Division  Naval  Surface  Warfare  Center  (CDNSWC),  Bethesda,  Maryland. 
T-ACS  Auxiliary  Crane  ships  are  ran  by  the  Military  Sealift  Command  Ready  Reserve 
Force.  These  self-sustaining,  rapidly-deployable  ships  support  military  sea  transportation 
needs  and  are  extremely  useful  in  ports  that  have  limited,  damaged  or  undeveloped  port 
facilities.  The  functions  of  T-ACS  ships  are  to  lift  and  transfer  various  loads  from  either 
themselves  or  adjacent  vessels  and  piers. 

A  total  of  251  different  model-scale  trial  runs  were  performed,  each  running 
approximately  8  minutes,  under  a  matrix  of  variable  conditions: 


Variable  Condition: 


Description: 


Ship  Configurations 


Config  I  T-ACS  in  Center 


Container  Ship  to  Port 


Lighter  Barges  to  Starboard 


Config  II  T-ACS  in  Center 


DDG-963  Class  Ship  to  Port 


Lighter  Barges  to  Starboard 


Ship  Heading 


Config  IE         T-ACS  in  Center 


Lighter  Barges  to  Port 


45  degree  increments  from  0  to  360  degrees  relative 


to  incoming  waves 


Sea  State 


Model-scale  sea  states  3,3  +  swell,  4  and  4  +  swell 


Boom  Slew  Angle 


Angles  of  0,  45,  90,  270  and  3 15  degrees 


Rider  Block  Location 


Various  positions  from  0  to  45.3  degrees 


Rider  Block  Inhaul  Angle       Various  angles  from  0  to  16  degrees 


Boom  Luff  Angle 


Angles  of  25,  29.6,  54.5  and  60  degrees 


A  composition  of  wave  height,  body  motion,  velocity  and  acceleration  data  for  the 
vessels  was  recorded  for  each  run  from  an  array  of  sensors  with  their  locations 
dependant   upon    the    ship   configurations    listed    above.  Results    from    the    213 

Configuration  I  trial  runs,  the  most  probable  ship  configuration,  were  provided  for 
frequency  response  analysis  of  the  raw  data  to  estimate  the  Full-Scale  Response 
Amplitude  Operators  (RAOs)  of  the  T-ACS  auxiliary  crane  ship.  An  RAO  is  basically  a 
measure  or  ratio  of  a  vessels  response  to  a  regular  wave  of  unit  amplitude  and  thus 
defined  accordingly: 


RAO  = 


amplitude  _of  _  response 
amplitude  _of  _  the  _  wave 


(Zubaly,1996,pp.  322) 


Tupper  expresses  the  method  of  studying  responses  in  a  seaway  in  the  following 
manner: 


"...  the  apparently  random  surface  of  the  sea  can  be  represented  by  the 
summation  of  a  large  number  of  regular  sinusoidal  waves,  each  with  its 
own  length,  height,  direction  and  phase. " 

(Tupper,  1996,pp.  104) 

"...  the  response  of  the  ship  in  such  a  sea  could  be  taken  as  the  summation 
of  its  responses  to  all  the  individual  wave  components.  Hence  the  basic 
building  block  for  the  general  study  of  motions  in  a  seaway  is  the  response 
to  a  regular  simusoidal  wave.  " 

(Tupper,  1996,pp.  104) 


These  concepts  strongly  identify  the  need  for  multiple  trial  runs  using  a  matrix  of  variable 
conditions.  The  study  of  responses  such  as  severe  rolling  in  a  beam  sea  or  excessive 
pitching  and  heaving  in  a  head  sea  provides  extremely  valuable  insight  into  a  vessels 
limitations.  Operating  or  load-handling  contraints  can  then  be  established  based  upon 
existing  as  well  as  predicted  conditions  to  prevent  material  damage,  downtime  and/or, 
even  worse,  human  injury. 
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II.       PROBLEM  FORMULATION 


A.        THEORY 


1.  Random  Process 

Determining  RAOs  begins  with  periodic  sampling  of  random  processes  such  as  a 
changing  waveheight  or  a  vessels  changing  responses  to  the  changing  waveheight  over  an 
interval  of  time.  In  this  case,  each  data  channel  was  sampled  at  a  rate  of  32.2  Hz  over  an 
interval  of  approximately  8  minutes  per  data  run.  This  sampling  produces  the  raw  data 
from  which  the  entire  RAO  determination  process  is  based  upon.  But,  before  the  raw 
data  can  be  used,  it  must  be  base-lined  to  remove  any  offset  or  bias  errors  inherent  in  the 
sampling  system  to  prevent  them  from  contaminating  the  process.  This  is  done  by 
determining  the  mean  of  the  data  points  for  a  given  channel  and  then  subtracting  this 
mean  from  each  data  point.  When  this  step  is  completed  for  each  data  channel,  the 
process  is  ready  to  continue. 


2.  Energy  Spectra  and  Power  Spectral  Densities  (PSDs) 

The  next  phase  in  determining  RAOs  is  to  generate  a  representative  energy 
spectrum  curve  for  each  data  channel  comprised  of  Power  Spectral  Densities  (PSDs)  vs. 
wave  frequencies  (omegas)  in  which  the  area  under  the  curve  represents  the  energy 
associated  with  the  changing  data.  (Papoulias,2000,pp.l09)  The  PSDs  for  this  curve  are 


derived  by  expansion  of  the  data  points  per  channel  into  a  Fourier  series  which 
approximates  the  shape  of  the  energy  spectrum  at  the  various  omegas. 
(Papoulias,2000,pp.  121)  The  mean  square  of  the  base-lined,  raw  data  record  within  a 
narrow  band  of  omega  centerd  at  omega  is  represented  as  follows: 


-2  - 

r  (CO)  =  S (CO) Act)         where  r  =  data _ mean _over_co_ band 


S(co)  =  energy  _  spectra 
Aco  =  omega  _band 
(Papoulias,2000,pp.  109) 


Thus,  as  an  integrity  check  along  the  way  which  is  utilized  later  in  the  PSD  and  RAO 
comparison  porgram,  the  mean  square  of  the  whole,  base-lined,  raw  data  record  should 
equal  the  area  under  the  entire  spectral  curve: 


/  (co)  =  jS(co)dco  (Papoulias,2000,pp.  109) 


3.  Response  Amplitude  Operators  (RAOs) 

Finally,  the  RAOs  are  determined  by  the  ratio  of  the  response  PSDs  over  the 
reference  waveheight  PSDs  as  follows: 


,DA-  ,2      PSDR 

\RAOJ   = -        where 

1       PSLL, 


RAOR  =  response  _  RAO 


B.         COMPARISON  TOOLS 


PSDR  =  response  _  PSD 


PSDW  =  waveheight  _  PSD 


(Papoulias,2000,pp.l24) 


1.  Spreadsheet  Results 

Spreadsheet  Results  are  selected  PSDs,  RAOs  and  Omegas  for  specific 
Configuration  I  trial  runs  that  were  initially  processed  by  Mr.  Dan  Hayden  of  CDNSWC 
using  a  spreadsheet  program.  Only  seven  data  channels  were  processed  per  run.  These 
channels  were  the  four  pitch  and  roll  channels  (channels  3-6)  as  well  as  the  three  relative 
TACS/Lighter  positions  (channels  28-30).  As  for  the  specific  data  runs,  only  19,  27,  57, 
67,  390  and  392  were  processed. 

The  Spreadsheet  Results  provided  a  reference  to  which  the  RAO  Program  results 
would  be  compared.  If  the  RAO  Program  results  were  consistant  with  each  of  the  seven 
processed  channels,  it  was  assumed  that  the  remaining  channels  would  be  correct  as  well. 
In  order  to  ensure  that  the  referenced  Spreadsheet  Results  were  correct,  three  additional 


MATLAB  methods  of  computing  PSDs  were  used  to  validate  both  the  Spreadsheet 
Results  and  the  RAO  Program. 


2.  MATLAB  Methods 

In  addition  to  the  FFT  method  in  the  RAO  Program,  three  additional  methods  of 
computing  PSDs  were  utilized  as  comparison  tools  to  ensure  validity  of  both  the 
Spreadsheet  Results  and  the  RAO  Program:  Welch,  Periodogram  and  Multi-taper 
Methods.  The  PSDs  are  computed  in  units  of  power  per  radians  while  the  frequencies  in 
Hz  are  later  converted  to  radians  per  second.  Although  there  are  many  similarities,  each 
method  performs  the  task  quite  differently. 

a.  Welch:  [PSDs,Freqs]  =  pwelch(x,nwindow,noverlap,nfftjs) 

The  selectable  options  when  utilizing  the  Welch  method  are  as  described 
as  follows: 


Option: Description 


x  Base-lined,  Raw,  Time-based  Data  Vector 

nwindow  Hamming   Window    Length   for   Modified   Periodograms 

(Default   of  []   divides  data  vector  into   8   equal-length 
windows  with  residual  data  points  discarded) 

noverlap  Number  of  Overlapping,  Windowed  Data  Points  (Default 

of  []  uses  50%  overlap  from  window  to  window) 
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nfft 


fs 


Length  of  the  Fourier  Transform  for  each  Window 


Data  Sampling  Frequency  (32.2  Hz) 


(SPT,welch) 


The  Welch  method  produces  PSDs  by  averaging  periodograms  of  overlapping, 
Hamming-windowed  sections  of  the  data  vector.  Hamming  windowing  reduces  possible 
sidelobes  in  the  spectral  estimate  in  order  to  reveal  the  presence  of  weaker  components  of 
the  signal  spectrum  that  may  otherwise  get  hidden.  (Marple,1987,pp.  132)  Specified- 
length,  discrete  Fourier  transforms  (DFTs)  of  the  overlapping  windows  are  computed  as 
follows: 


N-\ 


X  (k  + 1)  =  ]T  x(n  +  1)W„*"        where 


n=0 


WN=e 


—  „-/'(2*/A') 


N  =  length(x) 

X=DFT 

(SPT,welch) 


Indexing  of  (n+1)  and  (k+1)  is  used  since  MATLAB  vectors  run  from  1  to  N 
instead  of  0  to  N-l.  (SPT,fft)  A  modified  periodogram  for  each  windowed  segment  is 
then  computed: 


(\  N^, 


ModifiedPeriodogram  =  S(ejC0)  = 


1=1 


^>TI      I2 


vw;<=> 


(SPT,periodogram) 


The    individual    periodograms    are    averaged    to    produce    just    one    representative 
periodogram.  (SPT,periodogram)  Finally,  the  PSDs  are  computed: 


PSDs  = 


(  1  "l 


(SPT,periodogram) 


follows: 


Periodogram:  [PSDs,Freqs]  =  periodogram(xywindow,nfftjs) 


The  selectable  options  when  utilizing  the  Periodogram  method  are  as 


Option: 


Description 


x  Base-lined,  Raw,  Time-based  Data  Vector 

window  Window  Coefficients  for  Modified  Periodogram  of  the 

Input  Matrix  (Used  for  "m"x"n"  matricies.    Default  of  [] 
implies  single  column  data  vector) 
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nfft  Length  of  the  Fourier  Transform 

fs  Data  Sampling  Frequency  (32.2  Hz) 

(SPT,periodogram) 


The  Periodogram  method  produces  PSDs  in  the  same  general  manner  as  the 
Welch  method  previously  discussed.  A  modified  periodogram  is  generated  which  leads 
directly  to  computation  of  the  PSDs.  The  differences  between  the  two  are  that  no 
windowing/segmenting  of  the  data  vector  occurs  and  only  one  periodogram  is  generated. 
(SPT,periodogram)  Thus,  the  Welch  method  is  essentially  a  refinement  of  the 
Periodogram  method. 


c.  Multi-taper:    [PSDs,Freqs]  -pmtm(x,nwynfftjs) 

The  selectable  options  when  utilizing  the  the  third,  Multi-taper,  method 
are  as  follows: 


Option: Description 


x  Base-lined,  Raw,  Time-based  Data  Vector 

nw  Determines    Number    of    Discrete    Prolate     Spheroidal 

Sequences    (n  =  2mv-l)    used    as    Data   Tapers    for    of 
Estimation  of  PSDs  (Default  of  []  is  nw  =  4 ) 
11 


nfft  Length  of  the  Fourier  Transform 

fs  Data  Sampling  Frequency  (32.2  Hz) 

(SPT,pmtm) 


The  Multi-taper  method  is  by  far  the  most  complex  of  the  three  comparison 
methods  for  production  of  PSDs.  This  method  combines  linear  and  nonlinear  modified 
periodograms  to  estimate  the  PSDs  by  computing  each  periodogram  using  a  sequence  of 
orthogonal  tapers  or  windows  in  the  frequency  domain  as  specified  from  the  discrete 
prolate  spheriodal  sequences.  (SPT,pmtm) 
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III.     PROGRAM  DEVELOPMENT 


A.         GENERAL  DESCRIPTION 

The  RAO  program  contains  a  sequence  of  events  utilizing  standard  MATLAB 
functions.  The  validation  methods  to  be  discussed  later  are  dependant  upon  the 
MATLAB  Signals  Processing  Toolbox  from  which  the  commands  are  comprised  of 
standard  functions  within  the  source  code.  Each  significant  event  of  the  RAO  program  is 
addressed  below: 


1.  Initial  Set-up 

Before  processing  a  data  run,  both  the  raw  data  disc  and  the  output  storage  disc 
must  be  placed  in  their  respective  designated  drives.  Designation  of  these  drives  must 
conform  to  the  system  being  used  for  processing.  In  order  to  start  the  RAO  program,  the 
operator  enters  the  desired  run  number  to  be  processed  as  a  MATLAB  command  in  one 
of  the  following  two  formats: 


Command  Format: Applicable  DTMB  Runs: 

rao(##)  19-99 

rao(###)  101-587 
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Once  this  has  been  initiated,  the  program  designates  which  drive  the  raw,  time -based  data 
is  to  be  read  from  as  well  as  which  drive  the  generated  jpg  format  figures  and  composite 
matrices  for  the  specific  run  are  to  be  saved  to.  A  floating-point  format  with  5  digits  is 
placed  into  effect  and  the  following  constants  are  established: 


Constant: 


Definition: 


ptavg  =  50 


window  =  500 


freq  =  32.2 


Number  of  Sequential  PSDs  used  for  Smoothing 


Number  of  PSDs  within  a  Rational  Freq  Range 


lambda  =  24. 1 75  Actual  Model  Scale  (1:24.175) 


Raw  Data  Sampling  Frequency  in  Hz 


Selection  of  the  "ptavg"  constant  is  the  unavoidable  compromise  between  display 
resolution  and  accuracy  which  will  be  addressed  later.  As  for  the  "window"  constant,  it 
envelops  the  first  500  smoothed  PSDs  which  encompass  a  reasonable  range  of  omegas  to 
be  addressed  under  full-scale  conditions.  The  remaining  constant,  "freq",  is  the 
frequency  at  which  the  raw,  time -based  data  was  sampled. 


%     Initial      Set-up- 
function   rao(run) 


sread 
swri te 


char( 'H: \ ' ) ; 
chart 'D: \ ' 


format    short    e 
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ptavg  =  50; 

window  =  500; 

lambda  =  24.175; 

freq  =  32.2; 

2.         Establish  Proper  String  Configuration  for  the  Data  Run 

Based  upon  the  run  number  entered  by  the  operator,  proper  string  configurations 
for  the  run  as  well  as  the  folder  to  which  results  are  to  be  saved  are  established.  The 
program  then  differentiates  between  a  two  digit  and  a  three  digit  run  number  so  that  the 
run  string  evolves  as  "R###"  and  the  folder  as  "###".  Broad  applicability  of  designators 
such  as  these  are  utilized  throughout  the  RAO  program. 


%    Establish  Proper  String  Configuration   for  Data   Run- 
if  run   <=   99 

si  =     char('RO'); 

s2  =    char ( ' 0 ' ) ; 

else 

si  =     char('R'); 

end 

srun  =    num.2  str  (run 

filename  =    str cat  (si, srun) ; 

if  run   <=   99 

folder    =     s treat (s2 , srun) ; 
else 

folder    =     srun; 
end 


3.  Determine  Number  of  Data  Channels  to  be  Processed 

The  run  number  entered  enables  identification  of  the  number  of  data  channels  to 
be  processed  due  to  the  particular  data  channel  configuration  for  the  run  as  follows: 
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DTMB  runs: Configuration: Data  Channels: 

19-503  I  46 

504  -  533  II  38 

534  -  587  III  24 


Each  data  channel  is  individually  labeled  per  configuration  as  listed  in  Appendixes  A,  B 
andC. 


%  Determine  Number  of  Data   Channels    to  be   Processed 

if  run   <=   503 

channels      =     46; 
elseif  run   <=   533 

channels      =     38; 
else 

channels      =     24; 
end 


4.  Establish  File-path  and  Load  the  Raw,  Time-based  Data 

Given  the  run  number  and  drive  designation  from  which  the  data  is  to  be  read, 

"filepath"  is  established  as  a  string  in  order  to  execute  the  MATLAB  command  "load"  to 

load  the  raw,  time-based  data  matrix  for  processing.     Once  the  matrix  is  loaded,  its 

dimensions  are  identified  utilizing  the  MATLAB  command  "size."     The  number  of 

columns    representing    the    individual    data    channels    is    designated    as    the    variable 

"channels"  and  the  number  of  rows  representing  the  number  of  data  points  collected  is 

designated  as  the  variable  "n",  both  of  which  will  be  used  throughout  the  program. 
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%     Establish   Filepath   and  Load    the  Raw,     Time-based  Data 

filepath   =      [sread, filename] ; 

load ( filepath) 

szdata         =      size (data) ; 

channels    =      szdata  (2) ; 

n  =      szdata (1) ; 


5.  Compute  both  Model-Scale  and  Full-Scale  Omega  Ranges 

Before  going  any  further,  the  Model-Scale  omega  range  is  computed  based  upon 
the  sampling  frequency,  total  number  of  data  points  and  the  "window"  length  designated 
earlier.  Dividing  the  Model-Scale  range  by  the  square  root  of  the  model  scale  "lambda" 
produces  the  Full-Scale  omega  range.  This  is  supported  by  the  basic  assumption  for 
presentation  of  motion  data  that  "Natural  periods  of  motion  vary  as  the  square  root  of  the 
linear  dimension."  (Tupper,1996,pp.l06)  In  addition,  Tupper  identifies  that  "In  watching 
model  experiments  the  motion  always  seems  rather  'rapid'  because  of  the  way  period 
changes  .  Thus,  a  1/25  model  will  pitch  and  heave  in  a  period  only  a  fifth  of  the  full- 
scale  ship. 


%  Compute  both  Model-Scale   and  Full-Scale  Omega   Ranges 
momega  =    2  *pi  *freq/n  *  (0 :  window- 1 )  ; 

f omega  =    momega/ sqrt  (lambda)  ; 


6.  Full-Scale  Scalers,  PSD  Units  and  RAO  Units  Library 

There  are  four  individual  sets  of  Full-Scale  scalers,  PSD  units  and  RAO  units  in 
the  library  associated  with  the  type  of  data  collected  by  the  a  particular  channel.  The 
combinations  are  as  follows: 
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Data  Type: Full-Scale  Scaler:  PSD  Units:  RAO  Units: 

Wave  Height  1  in  dimensionless 

Distance  12/lambda  deg  deg/in 

Velocity  12/lambdaA(3/2)  deg/s  deg/in  s 

Acceleration  12/lambdaA(2)  deg/sA(2)  deg/in  sA(2) 


%     Full-Scale   Scalers,    PSD  Units   and  RAO  Units  Library- 

fsl  =  1 ; 

pul  =  '(in)'; 

rul  =  ' (dimensionless)  '  ; 

fs2  =  12/lambda; 

pu2  =  '  (deg)  '  ; 

ru2  =  ' (deg/in)  '  ; 

fs3  =  12/lambda"  (3/2)  ; 

pu3  =  '  (deg/s)  '; 

ru3  =  ' (deg/in  s)  ' ; 


fs4  =  12 /lambda" 2; 
pu4  =  '  (deg/s"2)  ' ; 
ru4      =      '(deg/in   s"2)  ' ; 


1.  Looping/Processing  of  Data  Channels 

It's  now  time  to  begin  looping  to  process  the  "channels"  sequentially.  From  the 
first  data  channel  to  the  last,  each  is  assigned  a  Full-Scale  Scaler,  "fscaler";  a  PSD  Units, 
"psdunits";  an  RAO  Units,  "raounits";  and  a  Channel  Title,  "chtitle".  These  assignments 
will  be  utilized  later  for  labeling  the  generated  jpg  figures. 
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%   Looping/Processing  of  Data   Channels- 
for   ch   =   1 : channels 

if  channels==4  6 

if  ch==l 

f  scaler        =  fsl; 

psdunits      =  pul; 

raounits      =  rul; 

chtitle        =  'Wave  Ht   Bow'; 
elseif  ch==2 

f  scaler        =  fsl; 

psduni ts      =  pul  ; 

raounits      =  rul; 

chtitle        =  ' Sonix  Sonic'; 


else   if  ch==46 

f  scaler        =  fs4; 

psdunits      =  pu4; 

raounits      =  ru4; 

chtitle        =  'Lghtrl-PBow  TvAcc ' ; 
else 
end 

elseif  channels==38 

if  ch==l 

f scaler        =  fsl; 

psdunits      =  pul; 

raounits      =  rul; 

chtitle        =  'Wave  Ht  Bow' ; 
elseif  ch==2 

f  scaler        =  fsl; 

psdunits      =  pul; 

raounits      =  rul; 

chtitle        =  'Sonix  Sonic'; 


elseif  ch==38 

f scaler        =    fs4; 
psdunits      =    pu4; 
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raounits  =  ru4 ; 

chtitle  =  'Lghtrl-PBowTvAcC  ; 
else 
end 

elseif  channel s==2 4 

if  ch==l 

f  scaler  =  fsl; 

psdunits  =  pul; 

raounits  =  rul ; 

chtitle  =  'Wave   Ht   Bow' ; 
elseif  ch==2 

f  scaler  =  fsl; 

psdunits  =  pul; 

raounits  =  rul; 

chtitle  =  ' Sonix  Sonic'; 


elseif  ch==24 

f  scaler  =     fs4; 

psdunits  =    pu4; 

raounits  =    ru4; 

chtitle  =      '  BoomTip-Vert  Ace'  ; 
else 
end 


end 


8.  Base-lining  the  Raw,  Time-based  Data 

The  actual  processing  of  data  begins  here  by  base-lining  or  removing  the  mean 
from  the  raw,  time-based  data.  Otherwise,  "Failure  to  remove  large  sample  means  or 
other  trends  in  the  data  may  result  in  distorted  or  biased  spectral  estimates." 
(Marple,1987,pp.l32) 
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%    Base-lining   the  Raw,    Time-based  Data- 
x  =      data  ( : , ch) ; 

x  =      x-mean  (x)  ; 


9.  FFT  Computation 

The  MATLAB  command  "fft"  returns  the  Fourier  transform  of  the  base-lined, 
column  data  matrix.  If  the  matrix  length  is  a  power  of  two,  the  "fft"  command  employs  a 
high-speed  fast  Fourier  transform  algorithm.  But,  this  is  not  always  the  case  since  the 
number  of  data  points  per  channel  varies  from  run  to  run  In  this  situation,  an  alternate 
mixed-radix  algorithm  finds  the  prime  factors  of  the  column  matrix  length  then  computes 
the  discrete  Fourier  transforms  of  the  shorter  sequences.  Either  way,  the  Fourier 
transform  of  the  base-lined  data  is  produced. 


%    FFT  Computation- 
Y  =     fft(x,n); 


10.        Model-Scale  PSD  Computation 

The  Model-Scale  PSD  is  computed  by  multiplying  the  elements  of  the  Fourier 
Transform  by  their  complex  conjugate  and  dividing  by  the  total  number  of  elements. 


%    Model-Scale   PSD  Computation- 
PSD      =     Y.  *conj  (Y)  /n; 
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11.        Moving-window  Averaging  of  Model-Scale  PSDs 

This  is  where  the  RAO  program  method  differs  significantly  from  the  tools  in 
MATLAB  previously  identified  by  offering  a  much  finer  control  over  the  resolution  vs. 
display  accuracy  issue.  An  envelope  initially  starts  with  half  the  "ptavg"  window  plus 
one  PSD  points,  averages  them  together  and  assigns  that  average  as  the  first,  smoothed 
PSD  point.  The  envelope  then  accepts  the  next  sequential  PSD  point,  averages  the  now 
larger  window  and  assigns  that  average  as  the  second,  smoothed  PSD  point.  The  process 
continues  until  the  envelope  reaches  it's  maximum  size  of  "ptavg"  plus  one  where  the 
latest  average  is  assigned  as  the  smoothed  PSD  point  for  the  middle  of  the  envelope. 
From  here  on  out,  the  envelope  accepts  the  next  sequential  PSD  point  while  dropping  the 
first  one,  thus  maintaining  a  fixed-length,  moving  window.  Refer  to  the  following 
example: 
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Example:         Let  "ptavg"  =  4 


[ Freshly  Produced  PSD  Data  Points I 

[#l       #2        #3        #4        #5        #6        #7        #8        #9        #101 


[  1st  avg  PSD ]  Window  starts  at  half  "ptavg"  plus  one  (i.e.  2  +  1) 

[       2n    avg  PSD       ]  Window  increments  in  size  by  one 

[  3rd  avg  PSD  ]  Max  capacity  of  ptavg  +  1  (i.e.  4  +  1) 

[  4    avg  PSD  ]  Window  begins  shift  to  the  right 

th 

[  5    avg  PSD  ]  Window  continues  to  shift 

[  6th  avg  PSD  ]  Window  still  shifting 


At  this  point  in  the  RAO  program,  the  smoothed  PSDs  are  designated  as  a  column 
of  data  in  a  composite  PSD  matrix  with  the  column  designation  correlating  directly  to  the 
data  channel  being  processed.  Thus,  the  first  channel  processed  becomes  the  first  column 
of  PSD  data  in  the  composite  PSD  matrix.  As  for  the  chosen  window  length  of  50  plus 
one,  tests  of  trial  runs  using  incrementally  larger  point  average  envelopes  yielded  a 
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compromise  of  "ptavg"  =  50  which  offered  the  best  display  resolution  while  maintaining 
accuracy  in  comparison  to  the  spreadsheet  results  provided  (see  figs.  1-3) 


%    Moving-window  Averaging  of  Model-Scale  PSDs    for 

Smoothing 
half ptavg  =    ptavg/ 2; 
sum  =     0  ; 

for  k   =    1 : ptavg 

sum      =     sum    +   PSD(k,  1)  ; 
if  k  >  half  ptavg 

psd (k-half ptavg,  ch)    =    sum/k; 
else 
end 
end 

pts  =    window  +  halfptavg; 

ptavgplusone =    ptavg  +   1; 

J  =    Is 

for  k       =    ptavgplusone : pts 

sum  =     sum   +   PSD(ptavg+j ,  1)    -    PSD(j,l); 

psd (k-half ptavg, ch)    =    sum/ptavg; 

3  =j+l; 

end 

meansqr  =  mean (x. "2 ) ; 

psdarea  =  trap  (momega,psd  ( : ,  ch)  )  ; 

psdscaler  =  meansqr /psdarea, ■ 

psd(:,ch)  =  psd( : , ch) *psdscaler 
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R067CH3  Model-Scale  TACS  Roll  PSDs  "ptavg"  Comparison 
ptavg=.26vs..S0vs.  76  vs.  100  vs  Spreadsheet  PSD  Results 


'3  ■  :  ■*■ 

'■  Omega  (rad/s)  • 


Figure  1. 


R067CK3  Model-Scale  TACS  Roll  RAOs  "ptavg,"  Comparison 
piavg  =  26  vs.  50  vs.  76  vs.  100  vs.  Spreadsheet  PSD  Results 


Figure  2. 
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R067CK3  Fqll-ScaleTACS  Roll  RAOs  *ptavg"' Comparison 
ptavg  =  26  vs.  90  *s.  76  vs.  1 00     , 


Figure  3. 


0;6     ,,;  0.8 

'Omega  (rad/s) 


12.        Computation  of  Model-Scale  RAOs 

As  previously  addressed  in  the  theory,  RAOs  are  obtained  by  taking  the 
individual  square  roots  of  the  desired  PSDs  divided  by  the  seaway  or  wave  height  PSDs. 
With  the  first  data  channel  for  each  configuration  allocated  as  the  reference  wave  height, 
the  first  channel  of  RAOs  processed  will  naturally  be  at  unity,  a  constant  dimensionless 
value  of  one. 


%     Computation   of  Model-Scale  RAOs 

X  =      psd(: ,ch) ./psd(: , 1) ; 

mrao ( : , ch)  =      sqrt (X) ; 
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13.        Computation  of  Full-Scale  RAOs 

Computation  of  the  Full-Scale  RAOs  is  simply  a  process  of  multiplying  the 
Model-Scale  RAOs  by  the  Full-Scale  Scaler.  This  previously  selected  scaler  is  based 
upon  the  nature  of  the  data  sampled. 


%  Computation   of  Full-Scale  RAOs 

frao ( : , ch)  =    mrao ( : ,  ch)  *f scaler; 


14.        Plot  Model-Scale  PSDs  and  Save  Figure  in  jpg  Format 

A  plot  of  the  smoothed,  Model-Scale  PSDs  vs.  their  associated  Model-Scale 
omegas  is  generated  (see  fig.  4)  then  saved  as  a  figure  in  jpg  format  in  a  folder  dedicated 
to  the  particular  run  and  titled  accordingly  with  the  run  number,  channel  number  and  a 
"p"  for  Model-Scale  PSDs  label. 


%    Plot  Model-Scale  PSDs   and  Save  Figure   in   jpg  Format- 

clf 

plot  (momega  (1,  : )  ,psd(: ,  ch) ) 

grid 

sch      =    num2str (ch) ; 

title ([ 'Model-Scale   PSDs    for  DTMB  Run    'srun,  '    / 

Channel    ' ,sch,  ' :     ' chtitle] ) 
xlabel ( 'Omega    (rad/s) ' ) 
ylabel(['PSD    'psdunits] ) 
saveas (gcf,  [swrite,  folder,  ' \ ' ,  filename,  'C',sch,  'p'], 

'jpg')  ; 
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Figure  4. 


15.        Plot  Model-Scale  RAOs  and  Save  Figure  in  jpg  Format 

A  plot  of  the  Model-Scale  RAOs  vs.  their  associated  Model-Scale  omegas  is 
generated  (see  fig.  5)  then  saved  as  a  figure  in  jpg  format  in  a  folder  dedicated  to  the 
particular  DTMB  run  and  titled  accordingly  with  the  run  number,  channel  number  and  an 
"m"  for  Model-Scale  RAOs  label. 


%     Plot   Model-Scale  RAOs   and   Save   Figure   in   jpg  Format- 

clf 

plot  (momega  (1,  : )  ,mrao  ( :  ,  ch)  ) 

grid 

title ([ 'Model-Scale   RAOs    for  DTMB  Run    ' srun,  '    / 

Channel    ',sch,':     'entitle]) 
xlabel ( 'Omega    (rad/s)') 
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ylabel ( [ 'RAO    'raounits] ) 

saveas (gcf, [swrite, folder,  '  \  ' ,  filename,  'C  ,sch,  'm'], 

'jpg') ; 


Figure  5. 
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16.        Plot  Full-Scale  RAOs  and  Save  Figure  in  jpg  Format 

A  plot  of  the  Full-Scale  RAOs  vs.  their  associated  Full-Scale  omegas  is  generated 
(see  fig.  6)  then  saved  as  a  figure  in  jpg  format  in  a  folder  dedicated  to  the  particular  run 
and  titled  accordingly  with  the  run  number,  channel  number  and  an  "f"  for  Full-Scale 
RAOs  label. 
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%     Plot   Full-Scale  RAOs   and  Save  Figure   in   jpg  Format— 

elf 

plot (f omega  (1,  : ) , frao  ( : ,  ch) ) 

grid 

title ([ 'Full-Scale  RAOs   for  DTMB  Run    ' srun,  '    / 

Channel     ',sch,  ':     ' chtitle] ) 
xlabel (' Omega    (rad/s)') 
ylabel ( [ 'RAO    'raounits] ) 
saveas (gef ,  [swrite,  folder,  '  \  ',  filename,  'C',sch,  'f'], 

jpg ' )  ; 


Processing  of  the  first  data  channel  is  now  complete.  The  RAO  program  then 
loops  back  to  event  7  listed  above,  Looping/Processing  of  Data  Channels,  until  all  data 
channels  for  the  specific  run,  as  determined  by  the  data  configuration,  have  been 
processed  to  generate  the  Model-Scale  PSDs,  Model-Scale  RAOs  and  Full-Scale  RAOs. 


06  O.f 

Omega  (rad/s) 


Figure  6. 


30 


17.        Save  Omega  Matrices  as  Last  Column  of  Composite  PSD  and  RAO 
Matrices 

For  the  purpose  of  convenient  analysis,  the  single-column  omega  matrices  are 
tacked  on  as  the  final  columns  of  their  respective  PSD  or  RAO  composite  matrices.  This 
basically  eliminates  the  necessity  of  the  jpg  formatted  figures  by  allowing  direct 
comparison  between  any  two  columns  of  a  composite  matrix  simply  by  using  a  plotting 
function.  Thus,  PSD  vs.  PSD,  PSD  vs.  omega,  RAO  vs.  RAO  and  RAO  vs.  omega  plots 
can  readily  be  produced 


%     Save  Omega  Matrices   as  Last   Column   of  Composite   PSD 

and  RAO  Matrices 
ch  =    ch   +   1; 

psd  ( : ,  ch)  =    momega  (1,  : )  '  ; 

mrao  ( : ,  ch)        =    momega  (1,  : ) 
frao ( : , ch)        =     f omega (1, : ) ' ; 


18.        Save  Composite  PSD  and  RAO  Matrices 

This  final  event  in  the  sequence  saves  the  composite  PSD  and  RAO  matrices 
along  with  the  jpg  formatted  figures  previously  saved  in  the  folder  dedicated  to  the 
particular  DTMB  run.  It's  now  time  to  process  another  DTMB  run. 


%     Save   Composite   PSD  and  RAO  Matrices 

save  (  [swrite,  folder,  '  \  ' ,  filename,  'psd.  txrt '],  'psd', 

' -ascii  ' ,  ' -  tabs ' ) ; 
save ( [swrite,  folder ,  ' \ ' ,  filename,  'mrao .  txt' ] ,  'mrao', 

' -ascii  ' ,  ' -  tabs ' ) ; 
save ( [swrite,  folder ,  ' \ ', filename,  ' frao .txt' ]  ,  'frao', 

' -ascii ' , ' -tabs ' ) ; 
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B.         VALIDATION 

Validation  of  the  RAO  Program  was  performed  by  utilizing  a  MATLAB 
comparison  program,  "psdandraocomparisons.m",  written  to  compare  the  various 
methods  of  resultant  PSDs;  Welch,  Periodogram,  Multi-taper,  Spreadsheet  and  RAO 
Program.  The  comparison  program  is  comprised  of  a  sequence  of  events  with  each 
explained  individually  as  applied  to  run  67: 


1.         Load  Applicable  Matrices 


Once  the  comparison  program  has  been  initiated,  the  following  matrices  are 


loaded: 


Matrix: 


Definition: 


R067 

Run67.txt 

R067psd.txt 

R067mrao.txt 

R067frao.txt 


Raw,  Time-based  Data  (All  46  Channels) 


Spreadsheet  Model-Scale  Omegas,  PSDs  and  RAOs 


RAO  Program  Model-Scale  PSDs  and  Omegas 


RAO  Program  Model-Scale  RAOs  and  Omegas 


RAO  Program  Full-Scale  RAOs  and  Omegas 
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%   Load  Applicable  Matrices- 
load  R067; 
load  Run67 . txt; 
load  R067psd.txt; 
load  R067mrao.txt; 
load  R067frao.txt; 


2.  Base-line  the  Raw  Time-Based  Data 

Prior  to  initiating  the  MATLAB  methods  of  Welch,  Periodogram  and  Multi-taper, 
the  raw,  time-based  data  for  both  Channel  1:  Wave  Ht  and  Channel  3:  TACS  Roll  are 
base-lined  for  continuity  just  as  it  is  done  in  the  RAO  Program. 


%  Base-line    the  Raw  Time-based  Data    from  R067 

x  =    data(:,l);  %   Wave  Ht   Data 

x  =    x   -  mean  (x)  ; 

z  =    data(:,3);  %   TACS  Roll   Data 

z  =     z    -  mean (z) ; 


3.  Compute  Wave  Ht  PSDs  using  MATLAB  Methods 

The  Welch,  Periodogram  and  Multi-taper  methods  are  used  to  compute  their 

Wave  Ht  PSDs  and  associated  Frequencies.     The  Frequencies  will  be  converted  to 

Omegas  later.    For  efficiency  of  computation,  a  maximum  of  8096  FFTs  is  designated. 

This  happens  to  be  the  largest  power  of  two  that  is  less  that  the  total  number  of  data 

points  per  channel  which  range  from  approximately   14,000  to   16,000.     Otherwise, 

MATLAB  uses  a  slower  approach  of  seeking  out  the  primes  of  the  total  number  of  data 

points  and  computing  the  FFTs  according  to  these  segments. 
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%   Computation   of 
[hwelch,  f welch] 
[hperio,  fperio] 
[hmul  ti ,  fmul  ti  ] 
meansqrx 
welchareah 
perioareah 
multiareah 
hwelchscaled 
hperioscaled 
hmultiscaled 


Wave  Ht   PSDs   and  Associated  Freqs- 

=  pwelch(x,  [],[]  ,8096,32.2)  ; 

=  periodogram (x, [] , 8096, 32 . 2) ; 

=  pmtm(x,  [],  8096,32.2) 

=  mean  (x.  *2 )  ; 

=  trap  (wwelch, hwelch) ; 

=  trap (wperio, hperio) ; 

=  trap (wmulti, hmulti) ; 

=  hwelch*meansqrx/ welchareah, ■ 

=  hperi o * meansqrx /peri oareah ; 

=  hmul ti *meansqrx/mul tiareah; 


4.         Compute  TACS  Roll  PSDs  using  MATLAB  Methods 


The  same  process  described  above  is  now  performed  for  the  TACS  Roll  PSDs. 


%   Computation   of 
[rwelch,  f welch] 
[rperio, fperio] 
[rmul  ti ,  fmul  ti  ] 
meansqrz 
welcharear 
perioarear 
multiarear 
rwelchscaled 
rperioscaled 
rmultiscaled 


TACS  Roll    PSDs  and  Associated  Freqs- 

=  pwelch(z,  [],[]  ,8096,32.2)  ; 

=  periodogram  (z,  [] , 8096, 32 .2) ; 

=  pmtm(z,  [],  8096,  32.2)  ; 

=  mean  (z . "2 ) ; 

=  trap  (wwel  ch ,  rwel  ch )  ; 

=  trap  (wperio,  rperio) ; 

=  trap (wmulti , rmul ti) ; 

=  rwelch*  meansqrz /welcharear ; 

=  rperio*meansqrz /perioarear ; 

=  rmul ti *meansqrz /mul tiarear ; 


5.  Convert  MATLAB  Method  Frequencies  to  Omegas 

Conversion  of  frequency  to  omega  is  done  simply  by  multiplying  the  frequency 
matrices  by  2tc. 
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%   Convert  Frequency   to   Omega 

wwel  ch  =    2  *pi  *  fwel  ch  ; 

wperi o  =2 *pi * fperi o ; 

wmulti  =    2*pi*fmulti; 


6.  Compute  Areas  under  the  PSD  Curves 

As  an  integrity  check  along  the  way,  the  mean  squares  of  the  base-lined  data  are 
compared  to  the  areas  under  the  PSD  curves.  At  this  point,  both  the  Spreadsheet  and 
RAO  Program  results  are  integrated  into  the  comparison  program  to  ensure  continuity 
between  the  various  methods  with  comparative  results  as  follows: 


a.  Wave  Ht  Comparison 


Method: 


Result: 


Mean  Square  of  Base-lined  Wave  Ht  Data 


Welch  Area  (rad/s) 


Periodogram  Area  (rad/s) 


Multi-taper  Area  (rad/s) 


Spreadsheet  Area  (rad/s) 


RAO  Program  Area  (rad/s) 


0.0124 
0.0124 
0.0124 
0.0124 
0.0126 
0.0124 
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From  the  results  listed  above  for  run  67,  the  Welch,  Periodogram,  Multi-taper  and 
RAO  Program  methods  agree  to  the  nearest  1/1 000th  when  rounded  to  5  decimal 
positions.  In  addition,  each  is  within   1 .6%  of  the  Spreadsheet  results. 


b.  TACS  Roll  Comparison 

Method: Result: 

Mean  Square  of  Base-lined  TACS  Roll  Data  0.006 1 

Welch  Area  (deg/s)  0.006 1 

Periodogram  Area  (deg/s)  0.0061 

Multi-taper  Area  (deg/s)  0.0061 

Spreadsheet  Area  (deg/s)  0.0059 

RAO  Program  Area  (deg/s)  0.006 1 


From  the  results  listed  above  for  run  67,  the  Welch,  Periodogram,  Multi-taper  and 
RAO  Program  methods  agree  to  the  nearest  1/1 000th  when  rounded  to  5  decimal 
positions.  In  addition,  each  is  within  3.3%  of  the  Spreadsheet  results. 
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%   Compute  Areas  under    the   Wave  Ht   PSD  Curves -- 

welchareah2  =     trap  (wwelch, hwelchscaled) ; 

perioareah.2  =     trap  (wperio,  hperioscaled)  ; 

multiareah.2  =     trap  (wmulti ,  hmultiscaled)  ; 


spreadareah  =     trap (run67 ( : , 2) , run67 ( : , 3) ) ; 

raoprogareah  =     trap  (R067psd( : ,  47) ,R067psd( : , 1) ) ; 

%   Compare  Areas   under    the   Wave  Ht   PSD  Curves   against 

meansqrs 
wavehtpsdareacomparison  =     [meansqrx  welchareah2 

perioareah2  multiareah2   spreadareah  raoprogareah] 

%   Compute  Areas   under    the   TACS  Roll   PSD  Curves 

welcharear2  =     trap  (wwelch,  rwelchscaled)  ; 

perioarear2  =     trap  (wperio ,  rperio scaled) ; 

multiarear2  =     trap  (wmulti, rmultiscaled) ; 

spreadarear  =     trap  (runS7 ( : , 2) ,  run67 ( : ,  4) ) ; 

raoprogarear        =     trap (R067psd( : , 47) ,  R067psd( : ,  3 ) ) ; 

%   Compare  Areas    under    the   TACS  Roll    PSD  Curves   against 

meqnsqrs 
rollpsdareacomparison  =     [meansgrz  welcharear2 

perioarear2  multiarear2   spreadarear  raoprogarear] 


7.         Graphical  Comparison  of  Model-Scale  PSDs 

Plots  of  Model-Scale  PSDs  from  all  five  methods  for  both  the  Wave  Ht  and 
TACS  Roll  are  generated  for  visual  comparison.  From  these  comparison  plots,  it  was 
determined  that  previously  discussed  "ptavg"  selection  of  50  was  a  very  good 
compromise  between  display  resolution  and  accuracy  when  compared  to  the  other  four 
methods  of  generating  PSDs.  (see  figs.  7-10) 
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%  Graphical  Model-Scale   Wave  Ht   PSDs   Comparison 

All    5  Methods 
figure (1) 

plot (wwelch,  hwelchscaled)  %   Welch  Method 

hold  on 

plot (wperio, hperioscaled)  %   Perio  Method 

plot  (wmulti, hmultiscaled)  %  Multi   Method 

plot (run67 ( : , 2 ) ,  run67 ( : ,  3 ) )  %   Spreadsheet 

plot (R067psd( : ,  47)  ,R067psd( : ,  1) ) ;  %  Rao   Program 

title ( 'R067CH01   Model-Scale   Wave  Ht   PSDs   Comparison' ) 
xlabel ( 'Omega    (rad/s) ' ) 
ylaheK'PSD    (dimesionless)  ' ) 
axis( [0    6   0    .015]) 
grid 


%  Graphical   Model-Scale   Wave  Ht   PSDs   Comparison 

Spreadsheet   vs.    RAO  Program   only 
figure (2) 

plot (run67 ( : , 2) , run67 ( : ,  3) )  %   Spreadsheet 

hold  on 

plot(R067psd( : ,47)  ,R067psd( :  ,1) )  %  Rao   Program 

title ('RO 67  CHOI   Model-Scale   Wave  Ht   PSDs   Comparison' ) 
xlabel ( 'Omega    (rad/s) ' ) 
ylabel ( 'PSD    (dimesionless) ' ) 
axis( [0    6   0    .  0065] ) 
grid 

%  Graphical  Model-Scale   TACS  Roll    PSDs   Comparison 

All    5  Methods 
figure (3) 

plot (wwelch, rwelchscaled)  %   Welch  Method 

hold  on 

plot (wperio, rperioscaled)  %   Perio  Method 

plot  (wmulti, rmultiscaled)  %  Multi   Method 

plot  (run67 ( : , 2)  ,  run67 ( : ,  4) )  %   Spreadsheet 

plot(R067psd( : ,47)  ,R067psd( :  ,3) ) ;  %   Rao   Program 

title ('R067CH03   Model -Scale   TACS  Roll    PSDs 

Comparison ' ) 
xlabel ( 'Omega    (rad/s) ' ) 
ylabel ( 'PSD    (deg)  ' ) 
axis(  [0    6   0    .02]) 
grid 
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%  Graphical    TACS  Roll    PSDs   Comparison- 
Spreadsheet   vs.    RAO  Program  only- 
figure  (4) 

plot (run67  ( : , 2) , run67 ( : ,  4) ) 
hold  on 

plot(R067psd(:  ,47)  ,R067psd(:  ,3)  ) 
title  ( 'R067CH03  Model-Scale   TACS  Roll   PSDs 

Comparison ' ) 
xlabel ( 'Omega    (rad/s) ' ) 
y label ('PSD    (deg)  ' ) 
axis(  [0    6   0    .009]) 
grid 


%   Spreadsheet 
%  Rao   Program 


0.015 


R067CH01  Model-Scats  Wave  Ht  PSOs  Comparison 
WeJoh  vs.  Periodogram  vs.  Multi-taper  vs.  Spreadsheet  vs.  RAO  Program 


0.01  - 


0.005 


Figure  7. 
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RC67CH01  Model-Scale  Wave  Ht  PSDs  Comparison 
'  Spreadsheet  Results  vs.  RAO  Program  Results 
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Figure  9. 
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_x  10 


R067CH03  ModePScaJe  TACS  Roll  PSOs  Comparison 
Spreadsheat  vs.  RAO  Program  Only 


Figure  10. 


8.         Compute  Model-Scale  RAOs 

Just  as  it  is  performed  in  the  RAO  Program,  Model-Scale  RAOs  are  computed  for 
the  Welch,  Periodogram  and  Multi-taper  methods  by  taking  the  square  of  the  scaled 
TACS  Roll  PSDs  divided  by  the  scaled  Wave  Ht  PSDs. 


%   Compute  Model-Scale   TACS  Roll   RAOs- 

Rl 

welchmrao 

R2 

periomrao 

R3 

multimrao 


=  rwelchscaled. /hwelchscaled; 

=  sqrt(Rl); 

=  rperioscaled. /hperio scaled; 

=  sqrt  (R2) ; 

=  rmultiscaled. /hmultiscaled; 

=  sqrt(R3); 
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9.         Graphical  Comparison  of  Model-Scale  RAOs 

Plots  of  Model-Scale  TACS  Roll  RAOs  from  all  five  methods  are  generated  for 
visual  comparison  and  they  provide  very  satisfactory  results,  (see  figs.  11-12) 


~  •  R067CH03  MSJeXScaleTACS RSS  RAOs "CwrifHsfsoh"  V ".TT^ZZ^ 
Welchvs.  Periodogrann  vs.  Multi-taper  vs.-  Spreadsheet  vs:  RAO  Program 


Note  that  the  Spreadsheet  and  RAO  Program 
RAOs  are  much  smoother  while  the  others'  are 
in  their  onginal.  unsfrioothed  state 


Figure  1 1 , 


42 


RCB7CH03  Model;Sca)eTACS  Roll  RAOs  Comparison' 
■   -  Spreadsheet  'vs.1  RAO  Pf  bgram  Only 


Figure  12. 


Omega  (rad/s) 


10.       Scale  to  Full-Scale  RAOs  and  Omegas 

In  order  to  compare  Full-Scale  RAOS  and  Omegas,  the  Model-Scale  RAOs  and 
Omegas  are  scaled  accordingly  just  as  they  are  done  in  the  RAO  Program. 


%   Scale    to  Full-Scale  RAOs  and  Omegas 

welchfrao  =    welchmrao*12/24 . 175; 

periofrao  =    periomrao*12/24 . 175 ; 

multifrao  =    multimrao*12/24 . 175; 


wwelchfull 
wperiofull 
wmultifull 


=  wwelch/sqrt (24.175)  ; 
=  wperio/sqrt (24.175)  ; 
=    wmulti/sqrt (24.175)  ; 
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11.        Graphical  Comparison  of  Full-Scale  RAOs 

Finally,  the  Full-Scale  RAOs  and  Omegas  for  the  Welch,  Periodogram,  Multi- 
taper  and  RAO  Program  methods  are  plotted  for  visual  comparison,  (see  figs.  13-14) 


%  Graphical   Full-Scale   TACS  Roll   RAOs   Comparison 

4  Methods  Not   including  Spreadsheet 
figure (7 ) 

plot  (wwelchfull ,  welch frao)  %   Welch 

hold  on 

plot (wperio full, perio frao)  %   Periodogram 

plot (wmulti full, multi frao)  %  Multi- taper 

plot(R067frao(:,47),R067frao(:,3))  %  RAO  Program 

title  ('RO 67  CEO  3   Full-Scale   TACS  Roll   RAOs   Comparison' ) 
xlabel (' Omega    (rad/s)') 
ylabel ( 'RAO    (deg/in)  ' ) 
axis([0   1.4    0   3]) 
grid 

%   Graphical   Full-Scale   TACS  Roll   RAOs   Presentation 

RAO   Program   only- 
figure  (8) 

plot  (RO 67 frao  ( : ,47) ,R067frao( :  ,3) )  %  Rao  Program 

title ('R067CH03   Full-Scale   TACS  Roll   RAOs 

Presentation ' ) 
xlabel (' Omega    (rad/s)') 
ylabel ( 'RAO    (deg/in)  ' ) 
axis([0   1.4    0   1.4]) 
grid 
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' R067CH03F«Scafe TACS'RoB RAOS Comparison    T K 
Welches.  Periodogram vs.Mutti-taper vs. RAO  Program \ 


0.6  0.8 

Omega  (rad/s) 


Figure  13. 


'■RD57C^3rFiffl^a&"TSCyRWRAOs"P^ehladon" 
:    RAO  Program  only    • 


0.2  0.4 


0.6  0.8 

Omega  (rad/s) 


Figure  14. 
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IV.      CONCLUSIONS  &  RECOMMENDATIONS 


A.         CONCLUSIONS 

Due  to  the  validation  comparisons  performed  on  run  67  above,  the  RAO  Program 
performed  quite  admirably  and  accurately.  The  RAO  Programs  ability  to  accurately 
produce  Model-Scale  PSDs,  Model-Scale  RAOs  and  ultimately  Full-Scale  RAOs  for 
each  channel,  along  with  their  associated  Omegas,  attests  to  the  programs  performance. 
In  addition,  a  short  looping  program,  "raoloop.m",  enabled  up  to  seven  data  runs  to  be 
evaluated  sequentially. 


%   Short   Program   Title:     "raoloop.m" 

a     =     [19   21    23   25   27   29   31];  %   DTMB  Run  Numbers 

for   z   =   1:7 

rao(a  (l,z)  )  ; 
end 


The  seven  run  looping  limit  was  imposed  due  to  local  hardware  configuration 
limitations.  With  each  data  run  producing  approximately  12.5  megabytes  of  stored 
figures  and  matrices,  a  100  megabyte  zip  disc  could  safely  store  up  to  seven  data  runs. 
With  a  read/write  CD  drive,  the  looping  capability  would  have  been  much  greater. 


Upon  completion  of  processing  the  Configuration  I  runs,  the  composite  matrices 
(Model-Scale  PSDs,  Model-Scale  RAOs  and  Full-Scale  RAOs)  for  each  run  were 
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assembled  all  on  one  compact  disc  as  a  convenient  reference  from  which  any  two 
columns  of  results  may  be  plotted  against  each  other. 


B.         RECOMMENDATIONS 

An  opportunity  for  a  logical  continuation  of  this  thesis  will  soon  be  available. 
Trial  runs  using  the  actual,  Full-Scale  ship,  the  USNS  Grand  Canyon  State  (ex-SS 
President  Polk)  (T-ACS  3),  were  conducted  September  9th- 16th  off  the  coast  of  Southern 
California  near  Camp  Pendleton.  Data  recorded  from  the  Full-Scale  trial  runs  could  be 
processed  in  a  manner  similar  to  the  RAO  Program  and  used  as  a  comparison  tool  to 
either  disprove  or  further  prove  the  validity  of  results  from  the  RAO  Program.  If  RAO 
results  from  both  the  at-sea  Full-Scale  trial  runs  and  the  RAO  Program  model-basin 
evaluations  are  in  agreement,  a  multitude  of  combinations  can  then  be  evaluated  to 
predict  ship  motions  in  an  irregular  sea  by  using  the  principle  of  linear  superposition. 
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APPENDIX  A:      TEST  CONFIGURATION  I 


CDNSWC 
Run  # 

Date 

Run 

Storage 

Directory 

Reference 
Zero  Run 

NRL 
Run 
# 

Test 

Conzig 

# 

Wave 
Head 
ing 

Sea 

State 

Boom 
Slew 
Angle 

Rider 
Block 
Location 

Rider 
Block 

Inhaul 
Angle 

Boom 
Luff 
Angle 

Comments 

019 

15- 

Jul- 

97 

TASCA2 

017,020 

001 

i 

0 

3 

0 

26.7 

0 

54.5 

021 

15- 

Jul- 
97 

TASCA2 

020,022 

002 

i 

0 

3s 

0 

26.7 

0 

54.5 

023 

15- 

Jul- 
97 

TASCA2 

022,024 

003 

i 

0 

4 

0 

26.7 

0 

54.5 

025 

15- 

Jul- 

97 

TASCA2 

024,026 

004 

i 

0 

4s 

0 

26.7 

0 

54.5 

027 

16- 

Jul- 

97 

TASCA2 

026,028 

005 

i 

90 

3 

0 

26.7 

0 

54.5 

Extreme 
pendulation 
perpendicular  to 
boom. 

029 

16- 
Jul- 

97 

TASCA2 

028,030 

006 

i 

90 

3s 

0 

26.7 

0 

54.5 

Load  hit  boom  & 
tangled  in  tag 
lines. 

031 

16- 

Jul- 

97 

TASCA2 

030,032 

007 

i 

90 

4 

0 

26.7 

0 

54.5 

Load  tangled  in 
tag  lines. 

033 

16- 

Jul- 
97 

TASCA2 

032,034 

008 

i 

0 

3 

0 

0 

0 

54.5 

24  file  format 

035 

16- 

Jul- 

97 

TASCA2 

034,036 

009 

i 

0 

3s 

0 

0 

0 

54.5 

24  file  format 

037 

16- 

Jul- 

97 

TASCA2 

036,038 

010 

i 

0 

4 

0 

0 

0 

54.5 

039 

16- 
Jul- 

97 

TASCA2 

038,040 

011 

i 

0 

4s 

0 

0 

0 

54.5 

041 

16- 

Jul- 

97 

TASCA2 

040,042 

012 

i 

90 

3 

0 

0 

0 

54.5 

044 

16- 

Jul- 

97 

TASCA2 

043,045 

014 

i 

0 

3 

0 

36.5 

3.5 

54.5 

046 

16- 

Jul- 
97 

TASCA2 

045,047 

015 

i 

0 

3s 

0 

36.5 

3.5 

54.5 

048 

16- 

Jul- 

97 

TASCA2 

047,049 

016 

i 

0 

4 

0 

36.5 

3.5 

54.5 

050 

16- 

Jul- 

97 

TASCA2 

049,051 

017 

i 

0 

4s 

0 

36.5 

3.5 

54.5 

057 

16- 

Jul- 

97 

TASCA2 

056,058 

018 

i 

90 

3 

0 

36.5 

3.5 

54.5 

059 

17- 

Jul- 

97 

TASCA2 

058,060 

019 

i 

90 

3s 

0 

36.5 

3.5 

54.5 
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17- 

061 

Jul- 
97 

TASCA2 

060,062 

020 

i 

90 

4 

0 

36.5 

3.5 

54.5 

17- 

063 

Jul- 
97 

TASCA2 

062,064 

021 

i 

90 

4s 

0 

36.5 

3.5 

54.5 

17- 

065 

Jul- 

97 

TASCA2 

064,066 

022 

i 

90 

2 

0 

36.5 

3.5 

54.5 

17- 

067 

Jul- 
97 

TASCA2 

066,068 

023 

i 

0 

2 

0 

36.5 

3.5 

54.5 

17- 

069 

Jul- 

97 

TASCA2 

068,070 

024 

i 

0 

2 

0 

26.7 

3.5 

54.5 

17- 

Repeated  as  Run 

081 

Jul- 
97 

TASCA2 

080,082 

l 

90 

2 

0 

26.7 

3.5 

54.5 

81  due  to  lost 
NRL  Run  025. 

17- 

073 

Jul- 

97 

TASCA2 

072,074 

026 

i 

90 

2 

0 

26.7 

12.3 

54.5 

17- 

075 

Jul- 
97 

TASCA2 

074,076 

027 

i 

90 

3 

0 

26.7 

12.3 

54.5 

17- 

077 

Jul- 

TASCA2 

076,078 

028 

i 

90 

3s 

0 

26.7 

12.3 

54.5 

97 

17- 

079 

Jul- 
97 

TASCA2 

078,080 

029 

i 

0 

3 

0 

26.7 

12.3 

54.5 

17- 

081 

Jul- 

97 

TASCA2 

080,082 

025 

i 

90 

2 

0 

26.7 

3.5 

54.5 

24  file  format 

17- 

083 

Jul- 
97 

TASCA2 

082,084 

030 

i 

90 

2 

0 

0 

0 

54.5 

17- 

085 

Jul- 
97 

TASCA2 

084,086 

031 

i 

0 

2 

0 

0 

0 

54.5 

17- 

087 

Jul- 
97 

TASCA2 

086,088 

032 

i 

0 

2 

270 

0 

0 

54.5 

24  file  format 

17- 

089 

Jul- 

97 

TASCA2 

088,090 

033 

i 

90 

2 

270  . 

0 

0 

54.5 

17- 

091 

Jul- 
97 

TASCA2 

090,092 

034 

i 

90 

3 

270 

0 

0 

54.5 

17- 

093 

Jul- 
97 

TASCA2 

092,094 

035 

i 

90 

3 

270 

0 

0 

54.5 

17- 

095 

Jul- 
97 

TASCA2 

094,096 

036 

i 

90 

3 

270 

0 

0 

54.5 

18- 

Lost  Barge  Rel 

098 

Jul- 
97 

TASCA2 

097,099 

037 

i 

0 

3 

270 

0 

0 

54.5 

Mot  Dyna  Site 
during  run. 

18- 

100 

Jul- 

97 

TASCA2 

099,101 

038 

i 

90 

2 

270 

26.7 

12.2 

54.5 

18- 

102 

Jul- 
97 

TASCA2 

101,103 

039 

i 

0 

2 

270 

26.7 

12.2 

54.5 

104 

18- 

TA3CA2 

103, 105 

04  0 

i 

90 

3 

270 

26.7 

12.2 

54.5 

50 


Jul- 

97 

18- 

106 

Jul- 

TASCA2 

105,107 

041 

i 

0 

3 

270 

26.7 

12.2 

54.5 

97 

18- 

108 

Jul- 

TASCA2 

107,109 

042 

i 

90 

3s 

270 

26.7 

12.2 

54.5 
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270 

36.5 

12 

54.5 

444 

14- 

Aug- 

97 

TASCA7 

443,445 

188 

I 

180 

4  + 

270 

36.5 

12 

54.5 

446 

14- 

Aug- 

TASCA7 

445,447 

189 

I 

180 

4s 

270 

36.5 

12 

54.5 

24  File  format 
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97 

14- 

448 

Aug- 
97 

TASCA7 

447,449 

190 

I 

180 

2 

270 

0 

0 

54.5 

14- 

450 

Aug- 

97 

TASCA7 

449,451 

191 

I 

180 

3 

270 

0 

0 

54.5 

14- 

452 

Aug- 

97 

TASCA7 

451,453 

192 

I 

180 

3s 

270 

0 

0 

54.5 

14- 

454 

Aug- 
97 

TASCA7 

453,455 

193 

I 

180 

4 

270 

0 

0 

54.5 

14- 

456 

Aug- 

97 

TASCA7 

455,457 

194 

I 

180 

4  + 

270 

0 

0 

54.5 

14- 

458 

Aug- 
97 

TASCA7 

457 

195 

I 

180 

4s 

270 

0 

0 

54.5 

15- 

460 

Aug- 
97 

TASCA7 

459,461 

196 

I 

180 

2 

315 

0 

0 

25 

15- 

462 

Aug- 

97 

TASCA7 

461,463 

197 

T 

270 

2 

315 

0 

0 

25 

15- 

464 

Aug- 

97 

TASCA7 

463,465 

198 

I 

180 

3 

315 

0 

0 

25 

15- 

466 

Aug- 
97 

TASCA7 

465,467 

199 

I 

270 

3 

315 

0 

0 

25 

15- 

468 

Aug- 

97 

TASCA7 

467,469 

200 

I 

180 

3s 

315 

0 

0 

25 

15- 

470 

Aug- 

97 

TASCA7 

469,471 

201 

I 

270 

3s 

315 

0 

0 

25 

15- 

472 

Aug- 
97 

TASCA7 

471,473 

202 

I 

180 

4 

315 

0 

0 

25 

15- 

474 

Aug- 
97 

TASCA7 

473,475 

203 

I 

270 

4 

315 

0 

0 

25 

NRL  Data  Lost. 

15- 

476 

Aug- 

97 

TASCA7 

475,477 

204 

I 

180 

4  + 

315 

0 

0 

25 

NRL  Data  Lost. 

15- 

478 

Aug- 

97 

TASCA7 

477,479 

205 

I 

180 

4s 

315 

0 

0 

25 

15- 

480 

Aug- 

97 

TASCA7 

479,481 

206 

I 

270 

4s 

315 

0 

0 

25 

15- 

482 

Aug- 

97 

TASCA7 

481,483 

207 

I 

270 

4  + 

315 

0 

0 

25 

15- 

484 

Aug- 

97 

TASCA7 

483,485 

208 

I 

180 

2 

315 

29.25 

12 

25 

15- 

486 

Aug- 
97 

TASCA7 

485,487 

209 

I 

270 

2 

315 

29.25 

12 

25 

15- 

488 

Aug- 

97 

TASCA7 

487,489 

210 

I 

180 

3 

315 

29.25 

12 

25 

NRL  Data  Lost. 
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490 

15- 

Aug- 

97 

TASCA7 

489,491 

211 

I 

270 

3 

315 

29.25 

12 

25 

NRL  Format  2 . 

492 

15- 

Aug- 

97 

TASCA7 

491,493 

212 

I 

180 

3s 

315 

29.25 

12 

25 

494 

15- 

Aug- 
97 

TASCA7 

493,495 

213 

I 

270 

3s 

315 

29.25 

12 

25 

496 

15- 

Aug- 

97 

TASCA7 

495,497 

214 

I 

180 

4 

315 

29.25 

12 

25 

NRL  Format  2 . 

498 

15- 

Aug- 
97 

TASCA7 

497,499 

215 

I 

270 

4 

315 

29.25 

12 

25 

NRL  Format  2 . 

500 

15- 

Aug- 

97 

TASCA7 

499,501 

216 

I 

180 

4s 

315 

29.25 

12 

25 

NRL  Format  2 . 

502 

15- 

Aug- 

97 

TASCA7 

501 

217 

I 

270 

4s 

315 

29.25 

12 

25 
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APPENDIX  B.      TEST  CONFIGURATION  II 


CDNSWC 
Run  # 

Date 

Run 

Storage 

Directory 

Reference 
Zero  Run 

NRL 
Run 
# 

Test 

Config 

# 

Have 
Head 
ing 

Sea 
State 

Boom 
Slew 
Angle 

Rider 
Block 
Location 

Rider 
Block 
Inhaul 
Angle 

Boom 
Luff 
Angle 

Comments 

504 

15- 

Aug- 
97 

TACSA7 

503,505 

218 

Ii 

180 

3 

315 

5 

12 

25 

506 

15- 

Aug- 

97 

TACSA7 

505,507 

219 

Ii 

270 

3 

315 

5 

12 

25 

508 

15- 
Aug- 

97 

TACSA7 

507,509 

220 

Ii 

270 

3s 

315 

5 

12 

25 

510 

15- 

Aug- 

97 

TACSA7 

509,511 

221 

Ii 

270 

4 

315 

5 

12 

25 

512 

15- 

Aug- 

97 

TACSA7 

511,513 

222 

Ii 

270 

4s 

315 

5 

12 

25 

515 

16- 

Aug- 

97 

TACSA7 

514,516 

223 

Ii 

90 

3 

270 

5 

12 

25 

NRL  Load  Data 
Not  Locked 

517 

16- 

Aug- 

97 

TACSA7 

516,518 

224; 

Ii 

90 

3s 

270 

5 

12 

25 

NRL  Format  2 . 

519 

16- 

Aug- 

97 

TACSA7 

518,520 

225 

Ii 

90 

4 

270 

5 

12 

25 

521 

16- 

Aug- 

97 

TACSA7 

520,522 

226 

Ii 

90 

4s 

270 

5 

12 

25 

523 

16- 

Aug- 

97 

TACSA7 

522,524 

223 

Ii 

90 

3 

270 

5 

12 

25 

Rerun  to  get 
good  load  data. 

526 

16- 

Aug- 

97 

TACSA7 

525 

NCD 

Ii 

90 

4 

Load  Control 
Configuration  1 

527 

16- 

Aug- 

97 

TASCA8 

525 

NCD 

Ii 

90 

4 

Load  Control 
Configuration  2 

528 

16- 

Aug- 

97 

TASCA8 

525 

NCD 

Ii 

90 

4 

Load  Control 
Configuration  2 

529 

16- 

Aug- 

97 

TASCA8 

525 

NCD 

Ii 

90 

4 

Transition 
Control 

Configuration  2 
to  3 

530 

16- 

Aug- 

97 

TASCA8 

525 

NCD 

Ii 

90 

4 

Load  Control 
Configuration  3 

531 

16- 

Aug- 

97 

TASCA8 

525 

NCD 

Ii 

90 

4 

Dynamic  Load 

Control 

Demonstration 

532 

16- 

Aug- 

97 

TASCA8 

525 

NCD 

Ii 

90 

4 

Load  Control 
Configuration  4 

533 

16- 
Aug- 

97 

TASCA8 

525 

NCD 

Ii 

90 

4 

Dynamic  Load 
Demo  No  Rider 
Block 
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APPENDIX  C:      TEST  CONFIGURATION  III 


CDNSWC 
Run  # 

Date 

Run 

Storage 

Directory 

Reference 
Zero  Run 

NRL 
Run 
# 

Test 
Conf ig 

# 

Wave 
Head 
ing 

Sea 
State 

Boom 
Slew 
Angle 

Rider 
Block 
Location 

Rider 
Block 

Inhaul 
Angle 

Boom 
Luff 
Angle 

Comments 

17- 

537 

Aug- 
97 

TASCA9 

536,538 

NCD 

Iii 

0 

2 

45 

9.5 

16 

60 

17- 

539 

Aug- 

97 

TASCA9 

538,540 

227 

iii 

90 

2 

45 

9.5 

16 

60 

17- 

541 

Aug- 
97 

TASCA9 

540,542 

228 

iii 

0 

3 

45 

9.5 

16 

60 

17- 

543 

Aug- 

97 

TASCA9 

542,544 

229 

iii 

90 

3 

45 

9.5 

16 

60 

17- 

545 

Aug- 
97 

TASCA9 

544,546 

230 

iii 

0 

3s 

45 

9.5 

16 

60 

17- 

547 

Aug- 
97 

TASCA9 

546,548 

231 

iii 

90 

3s 

45 

9.5 

16 

60 

17- 

549 

Aug- 
97 

TASCA9 

548,550 

232 

iii 

0 

4 

45 

9.5 

16 

60 

17- 

551 

Aug- 

97 

TASCA9 

550,552 

233 

iii 

90 

4 

45 

9.5 

16 

60 

17- 

553 

Aug- 
97 

TASCA9 

552,554 

234 

iii 

0 

4  + 

45 

9.5 

16 

60 

17- 

555 

Aug- 
97 

TASCA9 

554,556 

235 

iii 

0 

4s 

45 

9.5 

16 

60 

17- 

557 

Aug- 
97 

TASCA9 

556,558 

NCD 

iii 

90 

3 

45 

5.625 

15 

60 

• 

17- 

559 

Aug- 
97 

TASCA9 

558,560 

236 

iii 

0 

2 

45 

39.5 

0 

60 

17- 

561 

Aug- 
97 

TASCA9 

560,562 

237 

iii 

90 

2 

45 

39.5 

0 

60 

17- 

563 

Aug- 
97 

TASCA9 

562,564 

238 

iii 

0 

3 

45 

39.5 

0 

60 

17- 

565 

Aug- 
97 

TASCA9 

564,566 

239 

iii 

90 

3 

45 

39.5 

0 

60 

17- 

568 

Aug- 
97 

TASCA9 

567,569 

240 

iii 

0 

2 

90 

45.3 

13 

60 

17- 

570 

Aug- 
97 

TASCA9 

569,571 

241 

iii 

90 

2 

90 

45.3 

13 

60 

17- 

572 

Aug- 
97 

TASCA9 

571,573 

242 

iii 

0 

3 

90 

45.3 

13 

60 

574 

17- 

TASCA9 

573,575 

243 

iii 

90 

3 

90 

45.3 

13 

60 
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Aug- 

97 

576 

17- 
Aug- 

97 

TASCA9 

575,577 

244 

iii 

0 

3s 

90 

45.3 

13 

60 

578 

17- 
Aug- 

97 

TASCA9 

577,579 

245 

iii 

90 

3s 

90 

45.3 

13 

60 

580 

17- 
Aug- 

97 

TASCA9 

579,582 

246 

iii 

0 

4 

90 

45.3 

13 

60 

581 

17- 

Aug- 

97 

TASCA9 

579,582 

247 

iii 

0 

4  + 

90 

45.3 

13 

60 

583 

17- 

Aug- 

97 

TASCA9 

582 

248 

iii 

90 

3 

90 

45.3 

13 

60 

240secs  Static  + 
Load  Control 
Demo 

584 

17- 
Aug- 

97 

TASCA9 

582 

249 

iii 

0 

2 

90 

0 

0 

60 

240secs  Static  + 
Load  Control 
Demo 

585 

17- 
Aug- 

97 

TASCA9 

582 

250: 

iii 

90 

2 

90 

0 

0 

60 

240secs  Static  + 
Load  Control 
Demo 

586 

17- 
Aug- 

97 

TASCA9 

582 

251 

iii 

90 

3 

90 

0 

0 

60 

200secs  Static  + 
Load  Control 
Demo 
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APPENDIX  D:      MASTER  CHANNEL  LIST:  CONFIGURATION  I 


DTMB  Runs  1-503  and  NRL  Runs  1-217 

Channel  # 

Channel  Title 

Units 

Final  File    Title 

01 

Wave  Ht  Bow 

inch 

r###ch01 

02 

Sonix  Sonic 

inch 

r###ch02 

03 

TACS  Roll 

degree 

r###ch03 

04 

TACS  Pitch 

degree 

r###ch04 

05 

TACS  Roll  Rt 

deg/sec 

r###ch05 

06 

TACS  PitchRt 

deg/sec 

r###ch06 

07 

TACS  Yaw  Rt 

deg/sec 

r###ch07 

08 

TACS  CG  ZAcc 

G 

r###ch08 

09 

TACS  CG  YAcc 

G 

r###ch09 

10 

TACS  CG  XAcc 

G 

r###chl0 

11 

TACS  CT  ZAcc 

G 

r###chll 

12 

TACS  CT  YAcc 

G 

r###chl2 

13 

Cntr  Roll 

degree 

r###chl3 

14 

Cntr  Pitch 

degree 

r###chl4 

15 

Cntr  CG  ZAcc 

G 

r###chl5 

16 

Cntr  CG  YAcc 

G 

r###chl6 

17 

Cntr  CG  XAcc 

G 

r###chl7 

18 

TACS  Cntr  Rel  X 

inch 

r###chl8 

19 

TACS  Cntr  Rel  Y 

inch 

r###chl9 

20 

TACS  Cntr  Rel  Z 

inch 

r###ch20 

21 

Lhtr  Roll 

degree 

r###ch21 

22 

Lhtr  Pitch 

degree 

r###ch22 

23 

Lhtr  Roll  Rt 

deg/sec 

r###ch23 

24 

Lhtr  PitchRt 

deg/sec 

r###ch24 

25 

Lhtr  CG  ZAcc 

G 

r###ch25 

26 

Lhtr  CG  YAcc 

G 

r###ch26 

27 

Lhtr  CG  XAcc 

G 

r###ch27 

28 

TACS  Lhtr  Rel  X 

inch 

r###ch28 

29 

TACS  Lhtr  Rel  Y 

inch 

r###ch29 

30 

TACS  Lhtr  Rel  Z 

inch 

r###ch30 

31 

Load  Mo  Orthognl 

inch 

r###ch31 

32 

Load  Mo  Parallel 

inch 

r###ch32 

33 

BoomTip-HorizAcc 

G 

r###ch33 

34 

BoomTip-Vert  Ace 

G 

r###ch34 

35 

Lghtr4-Strn  VAcc 

G 

r###ch35 

36 

Lghtr4-Bow  VAcc 

G 

r###ch36 

37 

Lghtr3-Bow  VAcc 

G 

r###ch37 
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38 

Lghtr2-Bow  VAcc 

G 

r###ch38 

39 

Lghtrl-Bow  VAcc 

G 

r###ch39 

40 

Lghtr2-SMdShpVAc 

G 

r###ch40 

41 

Lghtr2-PMdShpVAc 

G 

r###ch41 

42 

Lghtr2-CLMdSpVAc 

G 

r###ch42 

43 

Lghtrl-CLMdSpVAc 

G 

r###ch43 

44 

Lghtr4-StrnLngAc 

G 

r###ch44 

45 

Lghtr4-PStrnTvAc 

G 

r###ch45 

46 

Lghtrl-PBowTvAcc 

G 

r###ch46 
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APPENDIX  E.      MASTER  CHANNEL  LIST:  CONFIGURATION  II 


DTMB  Runs  504-533  and  NRL  Runs  218-226 

Channel  # 

Channel  Title 

Units 

Final  File  Title 

01 

Wave  Ht  Bow 

inch 

r###ch01 

02 

Sonix  Sonic 

inch 

r###ch02 

03 

TACS  Roll 

degree 

r###ch03 

04 

TACS  Pitch 

degree 

r###ch04 

05 

TACS  Roll  Rt 

deg/sec 

r###ch05 

06 

TACS  PitchRt 

deg/sec 

r###ch06 

07 

TACS  Yaw  Rt 

deg/sec 

r###ch07 

08 

TACS  CG  ZAcc 

G 

r###ch08 

09 

TACS  CG  YAcc 

G 

r###ch09 

10 

TACS  CG  XAcc 

G 

r###chl0 

11 

TACS  CT  ZAcc 

G 

r###chll 

12 

TACS  CT  YAcc 

G 

r###chl2 

13 

Lhtr  Roll 

degree 

r###ch!3 

14 

Lhtr  Pitch 

degree 

r###chl4 

15 

Lntr  Roll  Rt 

deg/sec 

r###chl5 

16 

Lhtr  PitchRt 

deg/sec 

r###chl6 

17 

Lhtr  CG  ZAcc 

G 

r###chl7 

18 

Lhtr  CG  YAcc 

G 

r###chl8 

19 

Lhtr  CG  XAcc 

G 

r###chl9 

20 

TACS  LhtRelX 

inch 

r###ch20 

21 

TACS  LhtRelY 

inch 

r###ch21 

22 

TACS  LhtRelZ 

inch 

r###ch22 

23 

Load  Mo  Orthognl 

inch 

r###ch23 

24 

Load  Mo  Parallel 

inch 

r###ch24 

25 

BoomTip-HorizAcc 

G 

r###ch25 

26 

BoomTip-Vert  Ace 

G 

r###ch26 

27 

Lghtr4-Strn  VAcc 

G 

r###ch27 

28 

Lghtr4-Bow  VAcc 

G 

r###ch28 

29 

Lghtr3-Bow  VAcc 

G 

r###ch29 

30 

Lghtr2-Bow  VAcc 

G 

r###ch30 

31 

Lghtrl-Bow  VAcc 

G 

r###ch31 

32 

Lghtr2-SMdShpVAc 

G 

r###ch32 

33 

Lghtr2-PMdShpVAc 

G 

r###ch33 

34 

Lghtr2-CLMdSpVAc 

G 

r###ch34 

35 

Lghtrl-CLMdSpVAc 

G 

r###ch35 

36 

Lghtr4-StmLngAc 

G 

r###ch36 

37 

Lghtr4-PStrnTAcc 

G 

r###ch37 

38 

Lghtrl-PBowTvAcc 
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APPENDIX  F.      MASTER  CHANNEL  LIST:  CONFIGURATION  III 


DTMB  Runs  534-587  and  NRL  Runs  227-251 

Channel  # 

Channel  Title 

Units 

Final  File  Title 

01 

Wave  Ht  Bow 

inch 

r###ch01 

02 

Sonix  Sonic 

inch 

r###ch02 

03 

TACS  Roll 

degree 

r###ch03 

04 

TACS  Pitch 

degree 

r###ch04 

05 

TACS  Roll  Rt 

deg/sec 

r###ch05 

06 

TACS  PitchRt 

deg/sec 

r###ch06 

07 

TACS  Yaw  Rt 

deg/sec 

r###ch07 

08 

TACS  CG  ZAcc 

G 

r###ch08 

09 

TACS  CG  YAcc 

G 

r###ch09 

10 

TACS  CG  XAcc 

G 

r###chl0 

11 

TACS  CT  ZAcc 

G 

r###chll 

12 

TACS  CT  YAcc 

G 

r###chl2 

13 

DD963  Roll 

degree 

r###chl3 

14 

DD963  Pitch 

degree 

r###chl4 

15 

DD  963  CG  ZAcc 

G 

r###chl5 

16 

DD  963  CG  YAcc 

G 

r###chl6 

17 

DD  963  CG  XAcc 

G 

r###chl7 

18 

TACS!  DD963  RelX 

inch 

r###chl8 

19 

TACS1  DD963  RelY 

inch 

r###chl9 

20 

TACS1  DD963  RelZ 

Inch 

r###ch20 

21 

Load  Mo  Orthognl 

Inch 

r###ch21 

22 

Load  Mo  Parallel 

Inch 

r###ch22 

23 

BoomTip-HorizAcc 

G 

r###ch23 

24 

BoomTip-Vert  Ace 

G 

r###ch24 
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APPENDIX  G.      RAO  PROGRAM 


function   rao(run) 


%  Program  Title:  rao.m 

%  Referred  to  as  RAO  Program 

%  Must  enter  one  of  the  following  to  initiate  this  program: 

% 

%     rao(##)    for  DTMB  run  numbers  <100 

%     rao(###)     for  DTMB  run  numbers  >100 


%  Initial  Set-up 

sread  =  char ( ' H : \ ' ) ; 

swrite         =  char('D:\'); 


format  short  e 


ptavg 
window 
lambda 
f  req 


=  50; 

=  500; 

=  24.175; 

=  32.2; 


%  Establish  Proper  String  Configuration  for  Data  Run- 
if  run  <=  99 

si         =  char( 'R0' ) ; 

s2         =  char ( ' 0 ' ) ; 
else 

si        =  char ( 'R' ) ; 
end 


srun 
filename 


num2str (run)  ; 
strcat (si , srun) ; 


if  run  <=  99 
folder 
else 

folder 
end 


=  strcat (s2 , srun) ; 


=  srun; 


71 


%  Determine  Number  of  Data  Channels  to  be  Processed 

if  run  <=  5  03 

channels  =  46; 
elseif  run  <=  533 

channels  =  3  8; 
else 

channels   =  24; 
end 

%  Establish  Filepath  and  Load  the  Raw,  Time-based  Data- 

filepath       =   [sread, filename] ; 

load ( filepath) 

szdata  =      size (data); 

channels       =   szdata(2); 

n  =   szdata ( 1)  ; 

%  Compute  both  Model-Scale  and  Full-Scale  Omega  Ranges- 
momega  =  2*pi*freq/n* (0 : window-1) ; 

fomega         =  momega/ sqrt (lambda) ; 

%  Full-Scale  Scalers,  PSD  Units  and  RAO  Units  Library-- 

fsl  =  1; 

pul  =  ' (in) ' ; 

rul  =  ' (dimensionless) ' ; 

fs2  =  12 /lambda ; 

pu2  =  ' (deg) ' ; 

ru2  =  ' (deg/in) ' ; 

fs3  =  12 /lambda" (3/2) ; 

pu3  =  ' (deg/s) ' ; 

ru3  =  ' (deg/in  s) ' ; 

fs4  =  12/lambda"2; 

pu4  =  ' (deg/s"2) ' ; 

ru4  =   ' (deg/in  s"2) ' ; 

%  Looping/Processing  of  Data  Channels 

for  ch  =  1: channels 

if  channels==46 


if  ch==l 

fscaler  =  fsl 

psdunits  =  pul 

raounits  =  rul 
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chtitle    =  'Wave  Ht  Bow' ; 
elseif  ch==2 

f scaler    =  fsl; 


psdunits 

= 

pul; 

raounits 

= 

nil; 

chtitle 

= 

' Sonix  Sonic ' ; 

elseif  ch==3 

f scaler 

= 

fs2; 

psdunits 

= 

pu2; 

raounits 

= 

ru2; 

chtitle 

= 

'TACS 

Roll' ; 

elseif  ch==4 

f scaler 

= 

fs2; 

psdunits 

= 

pu2  ; 

raounits 

= 

ru2  ; 

chtitle 

= 

'TACS 

Pitch' ; 

elseif  ch==5 

f scaler 

= 

fs3; 

psdunits 

= 

pu3  ; 

raounits 

= 

ru3  ; 

chtitle 

= 

'TACS 

Roll  Rt' ; 

elseif  ch==6 

f scaler 

= 

fs3; 

psdunits 

= 

pu3  ; 

raounits 

= 

ru3  ; 

chtitle 

= 

'TACS 

Pitch  Rt' 

elseif  ch==7 

f scaler 

= 

fs3; 

psdunits 

= 

pu3  ; 

raounits 

= 

ru3  ; 

chtitle 

= 

'TACS 

Yaw  Rt ' ; 

elseif  ch==8 

f scaler 

= 

fs4; 

psdunits 

= 

pu4; 

raounits 

= 

ru4; 

chtitle 

= 

'TACS 

CG  ZAcc ' ; 

elseif  ch==9 

f scaler 

= 

fs4; 

psdunits 

= 

pu4; 

raounits 

= 

ru4; 

chtitle 

= 

'TACS 

CG  YAcc ' ; 

elseif  ch==10 

f scaler 

= 

fs4; 

psdunits 

= 

pu4; 

raounits 

= 

ru4; 

chtitle 

= 

'TACS 

CG  XAcc ' ; 
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elseif  ch==ll 

f scaler 

= 

fs4; 

psdunits 

= 

pu4; 

raounits 

= 

ru4; 

chtitle 

= 

'TACS 

CT  ZAcc ' ; 

elseif  ch==12 

f scaler 

= 

fs4; 

psdunits 

= 

pu4; 

raounits 

= 

ru4; 

chtitle 

= 

'TACS 

CF  YAcc ' ; 

elseif  ch==13 

f scaler 

= 

fs2; 

psdunits 

= 

pu2  ; 

raounits 

= 

ru2  ; 

chtitle 

= 

'Cntr 

Roll' ; 

elseif  ch==14 

f scaler 

= 

fs2; 

psdunits 

= 

pu2  ; 

raounits 

= 

ru2  ; 

chtitle 

= 

'Cntr 

Pitch' ; 

elseif  ch==15 

f scaler 

= 

fs4; 

psdunits 

= 

pu4; 

raounits 

= 

ru4; 

chtitle 

= 

'Cntr 

CG  ZAcc ' ; 

elseif  ch==16 

f scaler 

= 

fs4; 

psdunits 

= 

pu4; 

raounits 

= 

ru4; 

chtitle 

= 

'Cntr 

CG  YAcc ' ; 

elseif  ch==17 

f scaler 

= 

fs4; 

psdunits 

= 

pu4; 

raounits 

= 

ru4; 

chtitle 

= 

'Cntr 

CG  XAcc ' ; 

elseif  ch==18 

f scaler 

= 

fsl; 

psdunits 

= 

pul; 

raounits 

= 

rul; 

chtitle 

= 

'TACS 

Cntr  Rel  X' 

/ 

elseif  ch==19 

f scaler 

= 

fsl; 

psdunits 

= 

pul; 

raounits 

= 

rul  ; 

chtitle 

= 

'TACS 

Cntr  Rel  Y' 

/ 
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elseif  ch==20 

f scaler 

= 

fsl; 

psdunits 

= 

pul  ; 

raounits 

= 

rul  ; 

chtitle 

= 

'TACS 

Cntr  Rel  Z' 

/ 

elseif  ch==21 

f scaler 

= 

fs2; 

psdunits 

= 

pu2  ; 

raounits 

= 

ru2  ; 

chtitle 

= 

'Lhtr 

Roll' ; 

elseif  ch==22 

f scaler 

= 

fs2; 

psdunits 

= 

pu2; 

raounits 

= 

ru2; 

chtitle 

= 

'Lhtr 

Pitch' ; 

elseif  ch==23 

f scaler 

= 

fs3; 

psdunits 

= 

pu3  ; 

raounits 

= 

ru3  ; 

chtitle 

= 

'Lhtr 

Roll  Rt' ; 

elseif  ch==24 

f scaler 

= 

fs3; 

psdunits 

= 

pu3  ; 

raounits 

= 

ru3  ; 

chtitle 

= 

'Lhtr 

Pitch  Rt ' ; 

elseif  ch==25 

f scaler 

= 

fs4; 

psdunits 

= 

pu4; 

raounits 

= 

ru4; 

chtitle 

= 

'Lhtr 

CG  ZAcc ' ; 

elseif  ch==26 

f scaler 

= 

fs4; 

psdunits 

= 

pu4; 

raounits 

= 

ru4; 

chtitle 

= 

'Lhtr 

CG  YAcc ' ; 

elseif  ch==27 

f scaler 

= 

fs4; 

psdunits 

- 

pu4; 

raounits 

= 

ru4; 

chtitle 

= 

'Lhtr 

CG  XAcc ' ; 

elseif  ch==28 

f scaler 

= 

fsl; 

psdunits 

= 

pul; 

raounits 

= 

rul  ; 

chtitle 

= 

'TACS 

Lhtr  Rel  X' 

r 
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elseif  ch==29 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==30 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==31 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==32 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==33 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==34 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==35 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==36 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==37 
f scaler 
psdunits 
raounits 
chtitle 


fsl 
pul 
rul 
'TACS  Lhtr  Rel  Y' 

fsl 
pul 
rul 
'TACS  Lhtr  Rel  Z' 

fsl 
pul 
rul 
' Load  Mo  Orhognl ' ; 

fsl 
pul 
rul 
'Load  Mo  Parallel 

fs4 
pu4 
ru4 
'Boom  Tip-HorizAcc ' ; 

fs4 
pu4 
ru4 
' Boom  Tip-VertAcc ' ; 

fs4 
pu4 
ru4 
'Lghtr4-Strn  VAcc'; 

fs4 
pu4 
ru4 
'Lghtr4-Bow  VAcc'; 

fs4 

pu4 
ru4 
'Lghtr3-Bow  VAcc ' ; 
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elseif  ch==38 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==39 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==40 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==41 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==42 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==43 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==44 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==45 
f scaler 
psdunits 
raounits 
chtitle 


else 


end 


f scaler 
psdunits 
raounits 
chtitle 
%  End  of 


=  fs4 

=  pu4 

=  ru4 

=  'Lghtr2-Bow  VAcc ' ; 

=  fs4 

=  pu4 

=  ru4 

=  ' Lghtrl -Bow  VAcc' ; 

=  fs4 

=  pu4 

=  ru4 

=  ' Lghtr2-SMdShp  VAcc' ; 

=  fs4 

=  pu4 

=  ru4 

=  ' Lghtr2-PMdShp  VAcc ' ; 

=  fs4 

=  pu4 

=  ru4 

=  'Lghtr2-CLMdSp  VAcc' ; 

=  fs4 

=  pu4 

=  ru4 

=  ' Lghtrl -CLMdSp  VAcc' ; 

=  fs4 

=  pu4 

=  ru4 

=  'Lghtr4-StrnLng  Ace'; 

=  fs4 

=  pu4 

=  ru4 

=  'Lghtr4-PStrn  TvAcc ' ; 

=  fs4 

=  pu4 

=  ru4 

=  'Lghtrl-PBow  TvAcc ' ; 
Configuration  I  sub-Loop 


77 


elseif    channels==3 


if  ch==l 

f scaler 

= 

fsl; 

psdunits 

= 

pill; 

raounits 

= 

nil; 

chtitle 

= 

'Wave 

Ht  Bow' ; 

elseif  ch==2 

f scaler 

= 

fsl; 

psdunits 

= 

pul; 

raounits 

= 

rul  ; 

chtitle 

= 

' Sonix  Sonic ' ; 

elseif  ch==3 

f scaler 

= 

fs2; 

psdunits 

= 

pu2  ; 

raounits 

= 

ru2; 

chtitle 

= 

'TACS 

Roll' ; 

elseif  ch==4 

f scaler 

= 

fs2; 

psdunits 

= 

pu2  ; 

raounits 

= 

ru2  ; 

chtitle 

= 

'TACS 

Pitch' ; 

elseif  ch==5 

f scaler 

= 

fs3; 

psdunits 

= 

pu3  ; 

raounits 

= 

ru3  ; 

chtitle 

= 

'TACS 

Roll  Rt' ; 

elseif  ch==6 

f scaler 

= 

fs3; 

psdunits 

= 

pu3  ; 

raounits 

= 

ru3  ; 

chtitle 

= 

'TACS 

Pitch  Rt ' ; 

elseif  ch==7 

f scaler 

= 

fs3; 

psdunits 

= 

pu3  ; 

raounits 

= 

ru3  ; 

chtitle 

= 

'TACS 

Yaw  Rt ' ; 

elseif  ch==8 

f scaler 

= 

fs4; 

psdunits 

= 

pu4; 

raounits 

= 

ru4; 

chtitle 

= 

'TACS 

CG  ZAcc ' ; 
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elseif  ch==9 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==10 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==ll 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==12 
f sscaler 
psdunits 
raounits 
chtitle 

elseif  ch==13 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==14 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==15 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==16 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==17 
f scaler 
psdunits 
raounits 
chtitle 


fs4 
pu4 
ru4 
' TACS  CG  YAcc ' ; 

fs4 
pu4 
ru4 
' TACS  CG  XAcc ' ; 

fs4 
pu4 
ru4 
'TACS  CT  ZAcc' ; 

=  fs4; 

pu4; 

ru4; 

'TACS  CF  YAcc' ; 

fs2 
pu2 
ru2 
'Lhtr  Roll' ; 

fs2 
pu2 
ru2 
'Lhtr  Pitch' ; 

fs3 
pu3 
ru3 
'Lhtr  Roll  Rt' ; 

fs3 

pu3 
ru3 
'Lhtr  Pitch  RT' ; 

fs4 

pu4 
ru4 
' Lhtr  CG  ZAcc ' ; 


79 


elseif  ch==18 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==19 
fscaler 
psdunits 
raounits 
chtitle 

elseif  ch==20 
fscaler 
psdunits 
raounits 
chtitle 

elseif  ch==21 
fscaler 
psdunits 
raounits 
chtitle 

elseif  ch==22 
fscaler 
psdunits 
raounits 
chtitle 

elseif  ch==23 
fscaler 
psdunits 
raounits 
chtitle 

elseif  ch==24 
fscaler 
psdunits 
raounits 
chtitle 

elseif  ch==25 
fscaler 
psdunits 
raounits 
chtitle 

elseif  ch==26 
fscaler 
psdunits 
raounits 
chtitle 


fs4 
pu4 
ru4 
'Lhtr  CG  YAcC  ; 

fs4 
pu4 
ru4 
' Lhtr  CG  XAcc ' ; 

fsl 
pul 
rul 
'TACS  LhtRelX' ; 

fsl 
pul 
rul 
'TACS  LhtRelY' ; 

fsl 
pul 
rul 
'TACS  LhtRelZ' ; 

fsl 
pul 
rul 
' Load  Mo  Orthognl ' ; 

fsl 
pul 
rul 
'Load  Mo  Parallel'; 

fs4 
pu4 
ru4 
' BoomTip-HorizAcc ' ; 

fs4 
pu4 
ru4 
' BoomTip-Vert  Ace ' ; 


SO 


elseif  ch==27 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==28 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==29 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==30 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==31 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==32 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==33 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==34 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==35 
f scaler 
psdunits 
raounits 
chtitle 


fs4 
pu4 
ru4 
'Lghtr4-Strn  VAcc'; 

fs4 
pu4 
ru4 
'Lghtr4-Bow  VAcc'; 

fs4 
pu4 
ru4 
' Lghtr3-Bow  VAcc ' ; 

fs4 
pu4 
ru4 
'Lghtr2-Bow  VAcc'; 

fs4 
pu4 
ru4 
'Lghtrl-Bow  VAcc' ; 

fs4 
pu4 
ru4 
' Lghtr2-SMdShpVAc ' ; 

fs4 
pu4 
ru4 
' Lghtr2-PMdShpVAc ' ; 

fs4 
pu4 
ru4 
'Lghtr2-CLMdSpVAC ; 

fs4 
pu4 
ru4 
' Lghtrl-CLMdSpVAc ' ; 
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elseif  ch==36 
f scaler 
psdunits 
raounits 
chtitle 

elseif  ch==37 
f scaler 
psdunits 
raounits 
chtitle 

else 

f scaler 
psdunits 
raounits 
chtitle 


fs4 
pu4 
ru4 
' Lghtr4-StrnLngAc ' ; 

fs4 
pu4 
ru4 
'Lghtr4-PStrnTAcC  ; 

fs4 
pu4 
ru4 
'Lghtrl-PBowTvAcc'  ; 


end 


End  of  Configuration  II  sub-Loop 


else 


if  ch==l 

f scaler 

= 

fsl; 

psdunits 

= 

pul; 

raounits 

= 

rul; 

chtitle 

= 

'Wave 

Ht  Bow' ; 

elseif  ch==2 

f scaler 

=. 

fsl; 

psdunits 

- 

pul  ; 

raounits 

= 

rul; 

chtitle 

= 

' Sonix  Sonic ' ; 

elseif  ch= 

=  =  3 

f scaler 

= 

fs2; 

psdunits 

= 

pu2  ; 

raounits 

= 

ru2  ; 

chtitle 

= 

'TACS 

Roll' ; 

elseif  ch= 

=  =  4 

f scaler 

= 

fs2; 

psdunits 

= 

pu2  ; 

raounits 

- 

ru2  ; 

chtitle 

- 

'TACS 

Pitch' ; 

elseif  ch= 

=  =  5 

f scaler 

= 

fs3; 

psdunits 

= 

pu3  ; 

raounits 

= 

ru3  ; 

chtitle 

- 

'TACS 

Roll  Rt' 
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elseif  ch- 

==6 

f scaler 

=  fs3; 

psdunits 

=  pu3  ; 

raounits 

=  ru3  ; 

chtitle 

=   ' TACS 

Pitch  Rt 

elseif  ch= 

=  =  1 

f scaler 

=  fs3; 

psdunits 

=  pu3  ; 

raounits 

=  ru3  ; 

chtitle 

=   ' TACS 

Yaw  Rt ' ; 

elseif  ch= 

=  =  8 

f scaler 

=  fs4; 

psdunits 

=  pu4  ; 

raounits 

=  ru4  ; 

chtitle 

=   ' TACS 

CG  ZAcc' 

elseif  ch= 

--  =  9 

f scaler 

=   fs4; 

psdunits 

=  pu4  ; 

raounits 

=  ru4  ; 

chtitle 

=   ' TACS 

CG  YAcc' 

elseif  ch= 

==10 

f scaler 

=   fs4; 

psdunits 

=  pu4  ; 

raounits 

=  ru4  ; 

chtitle 

=   ' TACS 

CG  XAcC 

elseif  ch= 

:  =  11 

f scaler 

=   f  S4  ; 

psdunits 

=  pu4  ; 

raounits 

=  ru4  ; 

chtitle 

=   ' TACS 

CT  ZAcc' 

elseif  ch= 

==12 

f scaler 

=  fs4; 

psdunits 

=  pu4  ; 

raounits 

=  ru4  ; 

chtitle 

=   ' TACS 

CF  YAcc' 

elseif  ch= 

:  =  13 

f scaler 

=  fs2; 

psdunits 

=  pu2  ; 

raounits 

=  ru2  ; 

chtitle 

=   'DD  963  Roll' ; 

elseif  ch= 

:  =  14 

f scaler 

=   f  s2  ; 

psdunits 

=  pu2  ; 

raounits 

=  ru2  ; 

chtitle 

=   'DD  963  Pitch' 
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elseif  ch=- 
fscaler  = 
psdunits  : 
raounits  = 
chtitle  = 
elseif  ch== 
fscaler  = 
psdunits  = 
raounits  = 
chtitle  = 
elseif  ch=- 
fscaler  = 
psdunits  = 
raounits  = 
chtitle  = 
elseif  ch== 
fscaler  = 
psdunits  = 
raounits  = 
chtitle  = 
elseif  ch==19 

fscaler  = 
psdunits  = 
raounits  = 
chtitle  = 
elseif  ch== 
fscaler  = 
psdunits  = 
raounits  = 
chtitle  = 
elseif  ch== 
fscaler  = 
psdunits  = 
raounits  = 
chtitle  = 
elseif  ch== 
fscaler  = 
psdunits  : 
raounits  = 
chtitle  = 
elseif  ch== 
fscaler  = 
psdunits  : 
raounits  : 
chtitle    = 


:15 

:   fs4 

:   PU4 

:  ru4 

'DD  9  63  CG  ZAcc' ; 
■16 

■■  fs4 
i  pu4 
:  ru4 

'DD  963  CG  YAcc' ; 
17 

fs4 

pu4 

ru4 

'DD  963  CG  XAcc' ; 
18 

fsl 

pul 

rul 

'TACS1  DD9  63  RelX' ; 

fsl 

pul 

rul 

'TACS1  DD9  63  RelY' ; 
20 

fsl 

pul 

rul 

'TACS1  DD9  63  RelZ' ; 
21 

fsl 

pul 

rul 

' Load  Mo  Orthognl ' ; 
22 

fsl 

pul 

rul 

'Load  Mo  Parallel'; 
23 
:  fs4 

:   PU4 

ru4 

' BoomTip-HorizAcc ' ; 


S4 


else 

fscaler   =  fs4 
psdunits   =  pu4 
raounits   =  ru4 
chtitle    =  'BoomTip-Vert  Ace'; 
end   %  End  of  Configuration  III  sub-Loop 

end   %  End  of  the  Overall  Configuration  Loop 

%  Base-lining  the  Raw,  Time-based  Data 

x  =   data ( : , ch) ; 

x  =   x-mean  (x)  ; 


%  FFT  Computation 

Y  =  fft (x,n) ; 


%  Model-Scale  PSD  Computation 

PSD  =  Y.*conj (Y) /n; 

%  Moving-window  Averaging  of  Model-Scale  PSDs  for 

smoothing 

halfptavg       =  ptavg/2; 
sum  =  0 ; 

for  k  =  l:ptavg 

sum        =  sum  +  PSD(k,l); 
if  k  >  halfptavg 

psd(k-halfptavg, ch)  =  sum/k; 
else 
end 
end 

pts  =  window  +  halfptavg; 

ptavgplusone    =  ptavg  +  1; 

D  =  i; 

for  k   =  ptavgplusone : pts 

sum  =  sum  +  PSD(ptavg+j , 1)  -  PSD(j,l); 

psd (k-half ptavg, ch)  =  sum/ptavg; 

J  =  j  +  1; 

end 

meansqr  =  mean(x.A2)  ; 

psdarea  =  trap (momega, psd (:, ch)  )  ; 

psdscaler  =  meansqr /psdarea; 

psd(:,ch)  =  psd ( : , ch) *psdscaler; 


85 


%  Computation  of  Model-Scale  RAOs 

X  =   psd( : ,ch) . /psd( :  ,1)  ; 

mrao(:,ch)      =   sqrt (X) ; 

%  Computation  of  Full-Scale  RAOs 

frao(:,ch)      =  mrao (:, ch) *f scaler ; 

%  Plot  Model-Scale  PSDs  and  Save  Figure  in  jpg  Format 

elf 

plot (momega (1, : ) ,psd( : , ch) ) 

grid 

sch  =  num2str (ch) ; 

title ([ 'Model-Scale  PSDs  for  DTMB  Run  'srun,'  /  Channel 

' , sch, ' :  'entitle] ) 
xlabel (' Omega  (rad/s)') 
ylabel(['PSD  'psdunits]) 

saveas (gef , [swrite, folder, ' \ ' , filename, ' C , sch, 'p' ] , ' jpg 
' ) ;   %  D:\###\R###C##p. jpg 

%  Plot  Model-Scale  RAOs  and  Save  Figure  in  jpg  Format 

elf 

plot (momega ( 1 , : ) , mrao ( : , ch) ) 

grid 

title ([ 'Model-Scale  RAOs  for  DTMB  Run  'srun,'  /  Channel 

' ,sch, ' :  'entitle] ) 
xlabel ( 'Omega  (rad/s)') 
ylabel(['RAO  'raounits]) 

saveas (gef , [swrite, folder, ' \ ' , filename, ' C , sch, 'm' ] , ' jpg 
');   %  D:\###\R###C##m. jpg 

%  Plot  Full-Scale  RAOs  and  Save  Figure  in  jpg  Format 

elf 

plot ( f omega ( 1 , : ) , f rao ( : , ch) ) 

grid 

title ([ 'Full-Scale  RAOs  for  DTMB  Run  'srun,'  /  Channel 

' , sch, ' :  'entitle] ) 
xlabel ( 'Omega  (rad/s)') 
ylabel(['RAO  'raounits]) 

saveas (gef , [swrite, folder, ' \ ' , filename, ' C ' , sch, ' f ' ] , ' jpg 
' ) ;    %  D:\###\R###C##f . jpg 

end  %  End  of  all  Channel  Looping 
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%  Save  Omega  Matrices  as  Last  Column  of  Composite  PSD  and 

RAO  Matrices 

ch  =  ch  +  1; 


psd(:,ch)  =  momega ( 1 , 
mrao(:,ch)  =  momega ( 1 , 
frao(:,ch)        =  f omega (1, 


)  ' 
)  ' 


%  Save  Composite  PSD  and  RAO  Matrices 

%       D:\###\R###psd.txt 

%       D:\###\R###mrao.txt 

%       D:\###\R###frao.txt 

save ( [swrite, folder, ' \ ' , filename, 'psd. txt ' ] , 'psd' , ' -ascii 

' , '-tabs' ) ; 
save ( [swrite, folder, ' \ ' , filename, 'mrao . txt ' ] , 'mrao' , ' -ascii 

' , '-tabs' ) ; 
save ( [ swrite , folder,  ' \ ' , filename ,  ' f rao . txt ' ] ,  ' f rao ' ,  ' -ascii 

' , ' -tabs' ) ; 

%  End  of  RAO  Program 
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APPENDIX  H.      PTAVG  COMPARISON  PROGRAM 


%  Program  title:  ptavgcomparison.m 

clear 
clc 

%  Load  Spreadsheet  Results  for  run  67 

load  run67.txt; 

%  Load  RAO  Program  Results  for  ptavg  =  26,  50,  76  and  100-- 
load  R067psd26.txt, 
load  R067psd50.txt 
load  R067psd76.txt 
load  R067psdl00.txt; 

load  R067mrao26.txt 
load  R067mrao50.txt 
load  R0  67mrao7  6.txt 
load  R067mraol00.txt; 

load  R067frao26.txt 
load  R067frao50.txt 
load  R067frao76.txt 
load  R067fraol00.txt; 


%  Plot  Model-Scale  PSD  Comparison 

figured) 

plot (R067psd26 ( : , 47) ,R067psd26 ( : , 3) , . . . 

R067psd50 ( : , 47) ,R067psd50 ( : , 3) , . . . 

R0  67psd7  6( : ,47) ,R067psd7  6 ( : ,3) , . . . 

R067psdl00 ( : , 47) ,R067psdl00 ( : , 3) , . . . 
run67 ( : ,2) , run67 ( : ,4) ) 
title ( 'R067CH3  Model-Scale  TACS  Roll  PSDs  "ptavg" 

Comparison' ) 
xlabel ( 'Omega  (rad/s)') 
ylabel ( 'PSD  (deg) ' ) 
grid 
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%  Plot  Model-Scale  RAO  Comparison 

figure (2 ) 

plot (R067mrao26 ( 


R067mrao50 ( 
R067mrao76 ( 
R067mraol00 ( 
run67 ( : , 2 ) , run67 ( 


47) 

47) 
47) 
,47 
,5) 


R067mrao26 ( 
R067mrao50 ( 
R067mrao76 ( 
,R067mraol00 


,3)  ,  . 
,3)  ,  . 
,3)  ,  . 
(:,3) 


title ( 'R067CH3  Model-Scale  TACS  Roll  RAOs 

Comparison' ) 
xlabel ( 'Omega  (rad/s) ' ) 
ylabeK'RAO  (deg/in)  '  ) 
grid 


"ptavg" 


%  Plot  Full-Scale 
figure (3 ) 


RAO  Comparison- 


plot  (R067frao26  ( 
R067frao50( 
R067frao76( 


R067frao26 ( 
R067frao50 ( 
R067frao76( 


:  ,47) 
:,47) 
:,47) 

R0  67fraol00( : ,47)  , R067f raolOO ( :  ,3)  ) 
title ( 'R067CH3  Full-Scale  TACS  Roll  RAOs 

Comparison' ) 
xlabel (' Omega  (rad/s)') 
ylabeK'RAO  (deg/in)') 
grid 

%  End  of  ptavgcomparison  program 


"ptavg" 
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APPENDIX  I.       PSD  AND  RAO  COMPARISONS  PROGRAM 


%  Program  Title:  psdandraocomparisons .  m- 


%  Computational  and  Graphical  Comparison  between  the 
var ious  Me thods 


clear 
clc 

format  short 


Load  Applicable  Matrices' 


load  R0  67; 
load  run67.txt; 
load  R067psd.txt; 
load  R067mrao.txt; 
load  R0  67frao.txt; 


Raw,  Time-based  Data 
Spreadsheet  PSD  and  RAO  Results 
RAO  Program  Model-Scale  PSD  Results 
RAO  Program  Model-Scale  RAO  Results 
RAO  Program  Full-Scale  RAO  Results 


%  Base-line  the  Raw  Time-based-- 
x  =  data ( : , 1 ) ; 

x  =  x  -  mean (x)  ; 


%  Wave  Ht  Data 


data (  :  ,  3  )  ; 
z  -  mean ( z ) ; 


TACS  Roll  Data 


%  Computation  of 
[hwelch, fwelch] 
[hperio, fperio] 
[hmulti , fmulti ] 
meansqrx 
welchareah 
perioareah 
multiareah 
hwelchscaled 
hperioscaled 
hmulti scaled 


Wave  Ht  PSDs  and  Associated  Frequencies-- 

=  pwelch(x,  [],[]  ,8096,32.2)  ; 

=  periodogram(x,  [] , 8096, 32  .2)  ; 

=  pmtm(x,  [] ,8096,32.2)  ; 

=  mean (x. "2);  %  Wave  Ht  Data 

=  trap  (wwelch,  hwelch)  ;   %  Wave  Ht 

=  trap (wperio, hperio) 

=  trap (wmulti, hmulti) 

=  hwelch*meansqrx/ welchareah; 

=  hperio*meansqrx/perioareah; 

=  hmulti*meansqrx/multiareah; 
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%  Computation  of 
[rwelch, fwelch] 
[rperio, fperio] 
[rmulti, fmulti] 
meansqrz 
welcharear 
perioarear 
multiarear 
rwelchscaled 
rperioscaled 
rmultiscaled 


TACS  Roll  PSDs  and  Associated  Frequencies- 

=  pwelch(z, [],[], 8096,32.2) ; 

=  periodogram(z,  [] , 8096, 32  .2)  ; 

=  pmtm(z,  [] , 8096,32.2) ; 

=  mean (z . "2 ) ; 

=  trap (wwelch, rwelch) 

=  trap (wperio, rperio) 

=  trap (wmulti, rmulti) 

=  rwelch*meansqrz /welcharear, 

=  rperio*meansqrz /perioarear 

=  rmulti *meansqrz /multiarear 


TACS  Roll  Data 
TACS  Roll 


%  Convert  Frequency  to  Omega 

wwelch  =  2*pi*fwelch 

wperio  =  2*pi*fperio 

wmulti  =  2*pi*fmulti 


%  Compute  Areas  under  the  Wave  Ht  PSD  Curves- 

welchareah2 

perioareah2 

multiareah2 

spreadareah 

raoprogareah 


trap (wwelch, hwelchscaled) 
trap (wperio , hperioscaled) 
trap (wmulti , hmultiscaled) 
trap ( run67 ( : , 2 ) , run67 ( : , 3 ) ) 
trap(R067psd( : , 47) ,R067psd( 


1) 


%  Compare  Areas  under  the  Wave  Ht  PSD  Curves  against 

meansqrs 

wavehtpsdareacomparison     =  [meansqrx  welchareah2  . . . 
perioareah2  multiareah2  spreadareah  raoprogareah] 


Compute  Areas  under  the  TACS  Roll  PSD  Curves- 


welcharear2 
perioarear2 
multiarear2 
spreadarear 
raoprogarear 


trap (wwelch, rwelchscaled) 
trap (wperio , rperioscaled) 
trap (wmulti , rmultiscaled) 
trap(run67 ( : , 2) , run67 ( : , 4 ) ) ; 
trap(R067psd( : ,47) ,R067psd( : , 3 ) ) ; 


%  Compare  Areas  under  the  TACS  Roll  PSD  Curves  against 
meansqrs 

rollpsdareacomparison    =  [meansqrz  welcharear2  . . . 

perioarear2  multiarear2  spreadarear  raoprogarear] 
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%  Graphical  Model-Scale  Wave  Ht  PSDs  Comparison  of  All  5 
Me  thods 


figure (1) 

plot (wwelch, hwelchscaled) 

hold  on 

plot (wperio,hperioscaled) 

plot (wmulti , hmultiscaled) 

plot (run67 (: ,2) ,run67 ( : ,3) ) 

plot(R067psd( : ,47) , R067psd( 


1) 


%  Welch  Method 

%  Periodogram  Method 
%  Multi-taper  Method 
%  Spreadsheet  Results 
%  Rao  Program  Results 


title( 'R067CH01  Model-Scale  Wave  Ht  PSDs  Comparison') 

xlabel ( 'Omega  (rad/s)') 

ylabel('PSD  (dimesionless) ' ) 

axis( [0  6  0  .015] ) 

grid 


%  Graphical  Model-Scale  Wave  Ht  PSDs  Comparison 

%     Spreadsheet  vs .  RAO  Program  only- 
figure  (2 ) 

plot (run67 ( : , 2  )  , run67 ( : , 3 ) )  %  Spreadsheet  Results 

hold  on 
plot (R067psd( : , 47) ,R067psd( : , 1) )      %  Rao  Program  Results 

title( 'R067CH01  Model-Scale  Wave  Ht  PSDs  Comparison') 

xlabel (' Omega  (rad/s)') 

ylabel('PSD  (dimesionless)') 

axis( [0  6  0  .0065] ) 

grid 

%  Graphical  Model-Scale  TACS  Roll  PSDs  Comparison  of  All  5 

Me  thods 

figure (3 ) 

plot (wwelch, rwelchscaled) 

hold  on 

plot (wperio, rperioscaled) 

plot (wmulti , rmultiscaled) 

plot (run67 ( : ,2) ,run67 ( : ,4) ) 

plot (R067psd( : ,47) ,R067psd( : , 3 ) ) ; 


%  Welch  Method 

%  Periodogram  Method 
%  Multi-taper  Method 
%  Spreadsheet  Results 
%  Rao  Program  Results 


title ( 'R067CH03  Model-Scale  TACS  Roll  PSDs  Comparison') 

xlabel ( 'Omega  (rad/s)') 

ylabel ( 'PSD  (deg) ' ) 

axis (  [0  6  0  .02]  ) 

grid 
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%  Graphical  TACS  Roll  PSDs  Comparison 

%     Spreadsheet  vs  .  RAO  Program  only- 
figure  (4 ) 

plot (run67 ( : , 2) , run67 ( : , 4) )  %  Spreadsheet  Results 

hold  on 
plot (R067psd( : ,47) ,R067psd( : ,3) )      %  Rao  Program  Results 

title ( 'R067CH03  Model-Scale  TACS  Roll  PSDs  Comparison') 

xlabel (' Omega  (rad/s)') 

ylabeK  'PSD  (deg)  '  ) 

axis( [0  6  0  .009] ) 

grid 

%  Compute  Model-Scale  Roll  RAOs  from  MATLAB  Methods 

Rl  =  rwelchscaled. /hwelchscaled; 

welchmrao  =  sqrt(Rl); 

R2  =  rperioscaled. /hperioscaled; 

periomrao  =  sqrt(R2); 

R3  =  rmultiscaled. /hmultiscaled; 

multimrao  =  sqrt(R3); 

%  Graphical  Model-Scale  TACS  Roll  RAOs  Comparison  of  All  5 

Methods 

figure (5 ) 

plot (wwelch, welchmrao)  %  Welch  Method 

hold  on 

plot (wperio, periomrao)  %  Periodogram  Method 

plot (wmulti , multimrao)  %  Multi- taper  Method 

plot (run67 ( : , 2 ) , run67 ( : , 5 ) )  %  Spreadsheet  Results 

plot (R067mrao ( : , 47)  , R067mrao ( : , 3)  )  ;   %  Rao  Program  Results 

title( 'R067CH03  Model-Scale  TACS  Roll  RAOs  Comparison') 

xlabel ( 'Omega  (rad/s)') 

ylabeK 'RAO  (deg/in)') 

axis ( [0  6  0  10] ) 

grid 

%  Graphical  Model-Scale  TACS  Roll  RAOs  Comparison 

%     (Spreadsheet  vs.  RAO  Program  only) 

figure ( 6) 

plot (run67 ( : , 2) ,run67 ( : , 5) )  %  Spreadsheet  Results 

hold  on 

plot (R067mrao ( : , 47) , R067mrao ( : , 3) )     %  Rao  Program  Results 
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title ( 'R067CH03  Model-Scale  TACS  Roll  RAOs  Comparison') 

xlabel ( 'Omega  (rad/s)') 

ylabeK'RAO  (deg/in)  '  ) 

axis( [0  6  0  2.2] ) 

grid 


%  Scale  to  Full-Scale  RAOs  and  Omegas-- 
welchfrao  =  welchmrao*12/24 . 175 
periofrao  =  periomrao*12/24 . 175 
multifrao         =  multimrao*12/24 . 175 


wwelchfull 
wperiofull 
wmultifull 


=  wwelch/sqrt (24.175) 
=  wperio/sqrt (24 . 175) 
=  wmulti/sqrt (24 .175) 


Welch  Method 


%  Graphical  Full-Scale  TACS  Roll  RAOs  Comparison 

%     (4  Methods  Not  including  Spreadsheet) 

figure (7 ) 

plot (wwelchfull , welchf rao) 

hold  on 

plot (wperiofull , periofrao) 

plot (wmultifull, multifrao) 

plot (R067f rao ( : ,47) ,R067frao( :  ,3)  ) 


Periodogram  Method 
Multi-taper  Method 
RAO  Program  Results 


title ( 'R067CH03  Full-Scale  TACS  Roll  RAOs  Comparison') 

xlabel (' Omega  (rad/s)') 

ylabeK'RAO    (deg/in)') 

axis( [0    1.4    0    3]  ) 

grid 


%  Graphical  Full-Scale  TACS  Roll  RAOs  Presentation 

%     RAO  Program  only 

figure (8 ) 

plot (R067frao( : ,47) ,R067frao( : , 3) )     %  Rao  Program  Results 

title( 'R067CH03  Full-Scale  TACS  Roll  RAOs  Presentation') 

xlabel ( 'Omega  (rad/s)') 

ylabeK'RAO    (deg/in)') 

axis( [0    1.4    0    1.4] ) 

grid 

%  End  of  PSD  and  RAO  Comparison  Program 
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